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Colloidal dispersions are commonly encountered in everyday life and represent an important class 
of complex fluid. Of particular significance for many commercial products and industrial processes 
is the ability to control and manipulate the macroscopic flow response of a dispersion by tuning the 
microscopic interactions between the constituents. An important step towards attaining this goal 
is the development of robust theoretical methods for predicting from flrst-principles the rheology 
and nonequilibrium microstructure of well deflned model systems subject to external flow. In this 
review we give an overview of some promising theoretical approaches and the phenomena they seek 
to describe, focusing, for simplicity, on systems for which the colloidal particles interact via strongly 
repulsive, spherically symmetric interactions. In presenting the various theories, we will consider first 
low volume fraction systems, for which a number of exact results may be derived, before moving on 
to consider the intermediate and high volume fraction states which present both the most interesting 
physics and the most demanding technical challenges. In the high volume fraction regime particular 
emphasis will be given to the rheology of dynamically arrested states. 
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1. INTRODUCTION AND OVERVIEW 

Complex fluids exhibit a rich variety of flow behaviour 
which depends sensitively upon the thermodynamic con- 
trol parameters, details of the microscopic interparticle 
interactions and both the rate and specific geometry of 
] the flow under consideration. The highly nonlinear re- 
] spouse characteristic of complex fluids may be readily 
] observed in a number of familiar household products [if. 
For example, mayonnaise consists of a stabilized emul- 
sion of oil droplets suspended in water and behaves as 
a soft solid when stored on the shelf but flows like a 
i liquid, and is thus easy to spread, when subjected to 
] shear flow with a knife [^. This nonlinear viscoelastic 
J flow behaviour, known as shear-thinning, may be manip- 
1 ulated on the microscopic level by careful control of the 
3 oil droplet size distribution. In contrast, a dispersion of 
] corn-starch particles in water, at sufficiently high concen- 
] trations, exhibits a dramatic increase in shear viscosity 
] with increasing shear rate; a phenomenon called shear- 
I thickening [^, [j] . Even the familiar practical problem of 
] extracting tomato ketchup from a glass bottle presents 
a highly nonlinear flow. In this case the applied shear 
] stress, generally implemented by shaking, must exceed a 
i| critical value, the yield-stress, before the ketchup begins 
] to flow as desired. 
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Colloidal dispersions are a class of complex fluid which 
display all of the above mentioned nonlinear flow re- 
sponses 0. In addition to being of exceptional rele- 
vance for many technological processes, the considerable 
research interest in colloidal dispersions owes much to 
the existence of well characterized experimental systems 
for which the interparticle interactions may be tuned to 
relatively high precision (often possible by simply vary- 
ing the solvent conditions) Q . The ability to control the 
microscopic details of the colloidal interaction facilitates 
comparison of experimental results with theoretical cal- 
culations and computer simulations based on idealised 
models (see e.g. [3-Q)- In particular, the size of col- 
loidal particles makes possible light scattering, neutron 
scattering and microscopy experiments which provide in- 
formation inaccessible to experiments on atomic systems 
and which have enabled various aspects of liquid state 
theory to be tested in detail [l^l • 

The typical size of a colloidal particle lies in the range 
lOnm to Ifim and thus enables a fairly clear separation 
of length- and time-scales to be made between the col- 
loids and the molecules of the solvent in which they are 
dispersed. As a result, a reasonable first approximation 
is to represent the solvent as a continuum fluid, gener- 
ally taken to be Newtonian and thus characterised by 
a constant solvent viscosity (see FigH]). For suspended 
particles with a length-scale greater than approximately 
Ifim the continuum approximation of the solvent is com- 
pletely appropriate. However, this becomes questionable 
as the average size of the particles is reduced below a 
few nanometers, at which point the discrete nature of 
the solvent can no longer be ignored. Colloidal particles 
occupy an intermediate range of length-scales for which 
a continuum approximation for the solvent must be sup- 
plemented by the addition of first order Gaussian fluc- 
tuations (Brownian motion) about the average hydrody- 
namic fields describing the viscous flow of the continuum 
solvent. 

The Brownian motion resulting from solvent fluctua- 
tions not only plays an important role in determining the 
microscopic dynamics; it is essential for the existence of 
a unique equilibrium microstructure. With the impor- 
tant exception of arrested glasses and gels, the presence 
of a stochastic element to the particle motion allows a 
full exploration of the available phase space and thus 
enables application of Boltzmann-Gibbs statistical me- 
chanics to quiescent (and crgodic) colloidal dispersions. 
While the specific nature of the balance between Brown- 
ian motion, hydrodynamic and potential interactions de- 
pends upon both the observable under consideration and 
the range of system parameters under investigation, it 
is the simultaneous occurrence of these competing phys- 
ical mechanisms which gives rise to the rich and varied 
rheological behaviour of dispersions. Unfortunately, the 
complicated microscopic dynamics presented by disper- 
sions also serves to complicate the theoretical description 
of these systems |15|. 

The present review has been written with a number of 




FIG. 1: A schematic illustration of coarse graining as applied 
to colloidal dispersions. Continuum mechanics approaches 
treat the dispersion as a single continuum fluid (panel (a)), 
whereas a fully detailed picture is obtained by treating both 
colloids and solvent explicitly (panel (c)). The theoretical 
methods considered in this work operate at an intermediate 
level (panel (b)) in which the colloids are explicitly resolved 
but the solvent may be treated as a continuum. 



aims in mind. On one hand, we would like to present 
a relatively concise overview of the main phenomeno- 
logical features of the rheology of dispersions of spher- 
ical colloidal particles. In order to reduce the parameter 
space of the discussion, emphasis will be placed on the 
simple hard sphere model for which the space of con- 
trol parameters is restricted to two dimensions (volume 
fraction and flow rate). While both attractive colloids 
and the response to non-shear flows will be addressed, 
no attempt has been made to be comprehensive in this 
respect. Another primary aim of the present work is to 
provide an overview, within the context of the aforemen- 
tioned phenomenology, of microscopically motivated ap- 
proaches to the rheology and flow induced microstructure 
of colloidal dispersions. Although we will discuss some 
less well founded 'schematic model' approaches, the focus 
here is upon 'first-principles' theories which prescribe a 
route to go from a well defined microscopic dynamics to 
closed expressions for macroscopically measurable quan- 
tities. 

The formulation of a robust theory of dispersion rhe- 
ology from microscopic starting points constitutes a 
formidable problem in nonequilibrium statistical me- 
chanics. Although considerable progress has been made 
in this direction, a comprehensive constitutive theory 
analogous to that of Doi and Edwards for entangled linear 
polymers (23 - |25| remains to be found. At present there 
exist a number of alternative microscopic theoretical ap- 
proaches to dispersion rheology which, despite showing 
admirable success within limited ranges of the system 
parameters, have so far been unable to provide a uni- 
fied global picture of the microscopic mechanisms un- 
derlying the rheology of colloidal dispersions. Despite 
common starting points (the many-body Smoluchowski 
equation) the disparate nature of the subsequent approx- 
imations, each tailored to capture a particular physical 
aspect of the cooperative particle motion, make it diffi- 
cult to establish clear relations between different theoret- 
ical approaches. A goal of this work is thus to clarify the 
range of validity of the various theoretical approaches and 
to identify common ground. We note that the present 
work is well complemented by a number of recent reviews 
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addressing dispersion rheology from both experimental 
mil] and theoretical perspectives [MUllil- 

The paper is organized as follows: In section [2] we will 
discuss briefly some traditional continuum mechanics ap- 
proaches to rheology, both to give a feeling for the spirit 
of such work and to put into context some of the mi- 
croscopic results presented later (in section [7]) . In sec- 
tion [3] we will introduce and discuss in some detail the 
Smoluchowski equation defining the microscopic dynam- 
ics under consideration. In section 2] we will consider 
the equilibrium and non-equilibrium phase behaviour of 
hard-sphere colloids in the absence of flow, which is a 
necessary pre-requisite to the subsequent discussions. In 
section [5] we will give a brief overview of the relevant 
basic phenomenology of dispersion rheology, including 
the shear-thinning and shear-thickening of colloidal fluids 
and the yielding of colloidal glasses. In section [Hlwc will 
consider the various theoretical approaches to treating 
colloidal fluids under external flow. In particular, exact 
results for the microstructure and rheology of low vol- 
ume fraction systems and their (approximate) extension 
to finite fluid volume fractions are discussed in subsec- 
tions 16.21 and 16. 3[ respectively. In section [71 we consider 
the recently developed mode-coupling based approaches 
to the rheology of dense colloidal suspensions which en- 
able glass rheology to be addressed. Finally, in section [8] 
we will provide an outlook for future work and identify 
possible new avenues for theoretical investigation. 



2. CONTINUUM MECHANICS APPROACHES 

Rheology is primarily an experimental discipline. In- 
deed, one of the simplest experiments imaginable is to 
exert a force on a material in order to see how it deforms 
in response. More precisely, in a controlled rheological 
experiment one measures either the stress arising from a 
given strain or, more commonly, the strain accumulated 
following application of an applied stress. In practice, 
both stress controlled and strain controlled experiments 
are performed and provide complementary information 
regarding the response of a material sample. For the 
purpose of this review we will focus upon situations in 
which a homogeneous strain field is prescribed from the 
outset. The description of experiments for which macro- 
scopic stress is employed as a control parameter poses 
an enormous challenge for microscopically based theo- 
ries and demands careful consideration of the nontrivial 
mechanisms by which the applied stress propagates into 
the sample from the boundaries. 

Given the apparent complexity of any microscopic the- 
ory, it is quite natural to begin first at a more coarse- 
grained level of description in an effort to establish the 
general phenomenology and mathematical structure of 
the governing equations at the continuum level. Histor- 
ically, this methodology was pioneered by Maxwell in 
his 1863 work on viscoelasticity and continued to de- 
velop into the following century through the efforts of 



distinguished rheologists such as Rivlin and Oldroyd [30| . 
While much of this early work aimed to acheive a more 
fundamental mathematical understanding of viscoelastic 
response, strong additional motivation was provided by 
experiments on polymeric systems which exposed a large 
variety of interesting nonlinear rheological phenomena in 
need of theoretical explanation. Theoretical approaches 
to continuum rheology thus seek to obtain a constitutive 
equation relating the stress, a tensorial quantity describ- 
ing the forces acting on the system [3l|, to the deforma- 
tion history encoded in the strain tensor. 

The typical 'rational mechanics' approach to this prob- 
lem is to assume a sufficiently general integral or differen- 
tial constitutive relation between stress and strain and to 
then constrain this as much as possible via the imposition 
of certain exact or approximate macro scop ic symmetry, 
conservation and invariance principles [23l [30l |32| . The 
clear drawback to this methodology is that the entire 
particulate system is viewed as a single continuum field, 
thus losing any contact to the underlying colloidal inter- 
actions and microstructure ultimately responsible for the 
macroscopic response (see Fig. I). As a result, such con- 
stitutive theories are neither material-specific nor gen- 
uinely predictive in character. Despite these shortcom- 
ings, the continuum mechanics approach to rheology has 
attained a great level of refinement and can be applied to 
fit experimental data from a wide range of physical sys- 
tems d^llOl- Moreover, the experience gained through 
continuum mechanics modelling may well prove useful 
in guiding the construction of more sophisticated micro- 
scopic theories by providing constraints on the admiss- 
able mathematical form of the constitutive equations. 



2.1. The Lodge equation 

It is perhaps instructive to give an illustration of the 
spirit in which phenomenological constitutive relations 
may be constructed using continuum mechanics concepts. 
The example we choose is not only of intrinsic interest, 
but will also prove relevant to the discussion of a recent 
microscopically based theory of glass rheology [T7l - [l9| to 
be discussed in section[71 We consider a viscoelastic fluid 
subject to shear deformation with flow in the a;-direction 
and shear gradient in the y-direction (a convention we 
will continue to employ throughout the present work). 
Suppose that we wish to determine the infinitessimal 
shear stress dcxy at time t arising from a small strain 
increment d'y at an earlier time t' . As the material is 
viscoelastic, it is reasonable to assume that the influence 
of the strain increment d'^{t') = 'y{t')dt' on the stress at 
time t must be weighted by a decaying function of the 
intervening time t—t', in order to represent the influence 
of dissipative processes. Adopting a simple exponential 
form for the relaxation function it is thus intuitive to 
write 



da^y(t) = Goo exp 



t-t' 



iit')dt', 



(1) 
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where r is a relaxation time and Goo is an elastic constant 
(the infinite frequency shear modulus) . Assuming linear- 
ity, the total stress at time t may thus be constructed by 
summing up all of the infinitessimal contributions over 
the entire flow history, which we take to extend into the 
infinite past. We thus arrive at 



.it) 



dt' Go 



> exp 



t-t' 



Partial integration leads finally to 



O-xy(i) 



1 



dt' G{t-t')-f{t,t'), 



(2) 



(3) 



where G{t) — GooCxp[— t/r] is the shear modulus and 
'y{t,t') is the accumulated strain j{t,t') = J^,dsj{s). 
The simple integral relation ([3]) between shear stress and 
shear strain was first considered by Boltzmann. Indeed, 
the assumption that the stress increments ([1]) may be 
summed linearly to obtain the total stress is often re- 
ferred to as the 'Boltzmann superposition principle'. 

In order to extend ([3]) to a tensorial relation, i.e. a true 
constitutive equation, an appropriate tensorial general- 
ization of the accumulated strain j(t,t') must be iden- 
tified. For the spatially homogeneous deformations un- 
der consideration the translationally invariant deforma- 
tion gradient tensor E{t,t') transforms a vector ('mate- 
rial line') at time t' to a new vector at later time t via 
r{t) = E{t,t') •r(t'), where Eap = dra/drp. An alterna- 
tive nonlinear choice of strain measure is the symmetric 
Finger tensor B{t,t') = E{t,t')E'^ {t,t'). The Finger 
tensor contains information about the stretching of ma- 
terial lines during a deformation but is invariant with 
respect to solid body rotations of the material sample. 
For simple shear the Finger tensor is given explicitly by 



B = 




(4) 



where 7 = j{t,t'). The accumulated strain in the inte- 
grand of Eq. (jS]) can thus be identified as the xy element of 
B{t,t'). This suggests that the Boltmann integral form 
([3]) may be extended using the simple ansatz 



rrit) 



dt' B{t,t 



<t~t')/T 



(5) 



for the full stress tensor (see subsection 12.31 below for 
more justification of this nontrivial step). Equation ([5]) 
is known as the Lodge equation in the rheological liter- 
ature and is applicable in both the linear and nonlinear 
viscoelastic regime (23j . 



2.2. Upper convected Maxwell equation 

The assumption of an exponentially decaying shear 
modulus is generally attributed to Maxwell, who real- 
ized that this choice enabled an interpolation between 



a purely elastic response to deformations rapid on the 
time scale set by r and a viscous, dissipative response 
in the limit of slowly varying strain fields. In fact, the 
Lodge equation derived above is simply the integral form 
of a nonlinear (differential) Maxwell equation. In order 
to show this we first differentiate ([5]) to obtain 



Da- 

In 



1 Goo 

-cr = 1, 

T T 



(6) 



where we have introduced the uppcr-convected derivative 

^^^CT{t)-K{t)(T{t)-rT{t)K^{t), (7) 

and where the velocity gradient tensor K{t) is defined in 
terms of the deformation gradient tensor via 



^^E{t,t') ^ K{t)E{t,t') 



(8) 



For an incompressible material the stress is only deter- 
mined up to a constant isotropic term. Eq.Q may thus 
be expressed in an alternative form by first defining a 
new stress tensor 

S = <T-Gocl, (9) 
and substituting for cr in Eq.®. This yields 
nv 1 

^ + -S = Goo(/^(t)+K^(t)). (10) 
Dt T 

This differential form of the Lodge equation is known as 
the upper-convected Maxwell equation psj and is a non- 
linear generalization of Maxwell's original scalar model to 
the full deviatoric stress tensor. Historically, the uppcr- 
convected Maxwell equation was first proposed by 01- 
droyd HO] directly on the basis of Maxwell's differential 
form. 



2.3. Material objectivity 

The assumption that one can go from ([3]) to ([5]) on 
the basis of a single off-diagonal element appears at first 
glance to be rather ad hoc. On one hand, this choice can 
be justified retrospectively, using the fact that the Lodge 
equation ([5]) is derivable from a number of simple molec- 
ular models, e.g. the dumbell model for dilute polymer 
solutions [23|. However, from a continuum mechanics 
perspective ([S]) is the simplest generalization of which 
satisfies the 'principle of material objectivity'. This prin- 
ciple expresses the requirement that the constitutive re- 
lationship between stress and strain tensors should be 
invariant with respect to rotation of either the material 
body or the observer, thus preventing an unphysical de- 
pendence of the stress on the state of rotation. That this 
symmetry is an approximation becomes clear when con- 
sidering the material from a microscopic viewpoint: In 
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a noninertial rotating frame the apparent forces clearly 
lead to particle trajectories which depend upon the an- 
gular velocity. For many systems the neglect of these 
effects on the macroscopic response of the system is an 
extremely good approximation. For the overdamped col- 
loidal dynamics considered in this work inertia plays no 
role and the principle of material objectivity is exact [33j . 

Mathematically, it is straightforward to check whether 
or not a proposed tensorial constitutive equation is ma- 
terial objective. When subject to a time-dependent ro- 
tation R{t) the deformation gradient tensor transforms 
as 

E{t, t') = R{t)E{t, t')R^{t'), (11) 

where E is the deformation gradient in the rotating 
frame. The dependence of E upon the state of rota- 
tion arises because E contains information about both 
the stretching and rotation of material lines. Insertion of 
the transformed tensor (jlip into the constitutive equa- 
tion for the stress thus corresponds to a rotation of the 
material sample. Material objectivity is verified if the 
resulting stress tensor is given by 

&{t) ^ R{t)(T{t)R^{t). (12) 

As noted, the Finger tensor B contains only information 
about the stretching of material lines and transforms un- 
der rotation according to 

B{t, t') = R{t)B{t, t')R^{t). (13) 

The material objectivity of the Lodge equation ([S]), and 
thus the upper convected Maxwell equation (fTO|) . follows 
trivially from the fact that cr is a linear functional of B. 
Many phenomenological rheological models thus start by 
assuming a general functional dependence cr(t) = T[B] 
in order to guarantee a rotationally invariant theory. 

The vast majority of microscopically motivated theo- 
ries of dispersion rheology treat only a single scalar ele- 
ment of the stress tensor (generally the shear stress cTxy). 
Indeed, the rarity of microscopic tensorial constitutive 
theories may well be the primary reason for the apparent 
gap between continuum and statistical mechanical the- 
ories aiming to describe common phenomena. We will 
revisit the concept of material objectivity in section [7] 
when considering a recently proposed tensorial constitu- 
tive equation for dense dispersions. 

2.4. Beyond continuum mechanics 

In the last decade, significant progress has been made 
in understanding the response of colloidal dispersions to 
exernal flow on a level which goes beyond the fully coarse- 
grained phenomenological approaches of traditional con- 
tinuum rheology. Important steps towards a more re- 
fined picture have be en p rovided by studies based on 
mesoscopic models [26l l35l - [37| . However, while such phe- 
nomenological approaches can reveal generic features of 



the rheological response, they are not material specific 
and can therefore address neither the influence of the 
microscopic interactions on the macroscopic rheology nor 
the underlying microstructure, as encoded in the particle 
correlation functions. This deeper level of insight is pro- 
vided by fully microscopic approaches which start from a 
well defined particle dynamics and, via a sequence of ci- 
ther exact or clearly specified approximate steps, lead to 
closed expressions for macroscopically measureable quan- 
tities. The symmetry, invariance and conservation princi- 
ples used as input in the construction of continuum theo- 
ries, such as the material objectivity discussed in section 
12. 3i should then emerge directly as a consequence of the 
microscopic interactions. Such an undertaking clearly re- 
quires the machinery of statistical mechanics. 

Theories founded in statistical mechanics provide in- 
formation regarding the correlated motion of the con- 
stituent particles and arc thus capable, at least in princi- 
ple, of capturing non-trivial and potentially unexpected 
cooperative behaviour as exhibited by equilibrium and 
nonequilibrium phase transitions. This ability to cap- 
ture emergent phenomena is in clear contrast to contin- 
uum approaches where such physical mechanisms must 
be input by hand. An additional advantage of a sta- 
tistical mechanics based approach to rheology over the 
direct application of continuum mechanics is that impor- 
tant additional information is provided regarding the mi- 
crostructure of the system, as encoded in the correlation 
functions. It thus becomes possible to connect the con- 
stitutive relations to the underlying correlations between 
the colloidal particles and obtain microscopic insight into 
the macroscopic rheological response. Additional moti- 
vation to theoretically 'look inside' the flowing system is 
provided by developments in the direct visualization and 
tracking of particle motion in experiments on colloidal 
dispersions (confocal microscopy) [38l - l40{ . together with 
advances in the computer simulation of model systems 
under flow HH-lil. 

Although beyond the scope of the present work, we 
note that the influence of steady shear flow on glassy 
states has been addressed, albeit in an abstract setting, 
by generalized mean-field theories of spin glasses [13, S^j ■ 
Spin glass approaches have proved useful in describing 
the dynamical behaviour of quiescent systems [46| . In 
order to mimic the effect of shear flow a nonconservative 
force is introduced to bias the dynamics and break the 
condition of detailed balance characterizing the equilib- 
rium state [13] . While the abstract nature of these treat- 
ments certainly lends them a powerful generality, the lack 
of material specificity makes difficult a direct connection 
to experiment. 



3. MICROSCOPIC DYNAMICS 

Before addressing the phenomenology (section [5]) and 
approximate theories (sections [5] and [7]) of colloid rheol- 
ogy it is rewarding to first consider in detail the micro- 



scopic equation of motion determining the overdampcd 
colloidal dynamics. By a careful assessment of the funda- 
mental equation of motion a number of general observa- 
tions and comments can be made regarding the character 
of nonequilibrium states, solutions in special limits, im- 
portant dimensionless parameters and influence of hydro- 
dynamics, which arc independent of the specific system 
or approximation scheme under consideration. 

We consider a system consisting of N Brownian col- 
loidal particles interacting via spherically symmetric 
pairwise additive interactions and homogeneously dis- 
persed in an incompressible Newtonian fiuid of given 
viscosity. The probability distribution of the TV-particle 
configuration is denoted by ^{t) and satisfies the Smolu- 
chowski equation 



dt 



(14) 



where the probability flux of particle i is given by 



j. = v,(<)vl/(i) - ^ A, • [d, - /3F,)vI/(i), (15) 




2DqR t 



FIG. 2: The mean-squared-displacement of non-interacting 
colloidal particles in flow (x), gradient (y) and vorticity (z) 
directions as a function of time. The MSD in flow direction ex- 
liibits enhanced diffusion ('Taylor dispersion') for values of the 
shear strain greater than unity. Also shown are contour plots 
of the (non-normalised) probability distribution Pi{r,t)N{t) 
(see Ea ll6|l in the z = plane at times 2DoR^^t = 0.15, 1 and 
5, demonstrating shear induced anisotropy for 7 > 1 related 
to the onset of Taylor dispersion. 



where /3 = l/fc^T is the inverse temperature. The hy- 
drodynamic velocity of particle i due to the applied flow 
is denoted by Vi{t) and the diffusion tensor Dij de- 
scribes (via the mobility tensor Tij = fiDij) the hydro- 
dynamic mobility of particle i resulting from a force on 
particle j. The hydrodynamic velocity can be decom- 
posed into affine and particle induced fluctuation terms 
v,(i) =K{t)- r; vf (i), where vf (i) can be expressed in 
terms of the third rank hydrodynamic resistance tensor 
psj . The force on particle j is generated from the to- 
tal potential energy according to Fj ~ —djUN, where, in 
the absence of external fields, Un depends solely on the 
relative particle positions. The three terms contributing 
to the flux thus represent the competing effects of (from 
left to right in (jlSp ) external flow, diffusion and interpar- 
ticle interactions. 

While the Smoluchowski equation ([Til) is widely ac- 
cepted as an appropriate starting point for the treat- 
ment of colloidal dynamics, alternative approaches based 
on the Fokker-Planck equation have also been investi- 
gated [4^ . On the Fokker-Planck level of description the 
distribution function retains a dependence on the parti- 
cle momenta. Although this makes possible the treat- 
ment of systems with a temperature gradient (leading 
to thermophoretic effects), considerable complications 
arise when attempting to treat hydrodynamic interac- 
tions which make preferable the Smoluchowski equation. 

For the special case of monodisperse hard-spheres at 
flnite volume fractions under steady flow Ea. ([H)) can be 
numerically integrated over the entire fluid range using 
computationally intensive Stokesian dynamics simulation 
IS E^: mi . This simulation technique includes the full 
solvent hydrodynamics and provides a useful benchmark 
for theoretical approaches (for an overview of the com- 
puter simulation of viscous dispersions we refer the reader 



to [53 and references therein). While Stokesian dynam- 
ics simulations have focused primarily on simple shear, 
results have also been reported for extensional flow ge- 
ometries [5311. 



3.1. Non-interacting particles 

For the special case of non-interacting particles (Fj — 
0) equations ([H)) and describe the conflgurational 
probability distribution of an ideal gas under externally 
applied flow and may be solved analytically using the 
method of characteristics [s^ . For non- interacting par- 
ticles under steady shear the many-particle distribution 
function is given by a product of single particle functions 
^'({ri},t) = Pi(ri,t) X • • • X Pi(rN,t), where Pi is given 



by 



Nit) 



exp 



2/^(1 + ^ 



-jxy 



Nit) 



ADot 



(16) 



when the initial condition Pi (r, 0) = (5(0) is employed. 
The Normalization is given by N{t) = i4nDot)^il + 
(7t)^/12) and the strain by 7 = jt. Given a suitably 
localized initial density distribution Ea. (jl6p essentially 
describes the dispersion of a colloidal droplet in a solvent 
(e.g. ink in water) under shear, as is apparent from Fig|21 
which shows contour plots of the probability distribution 
at three different times for a given shear rate. 

Although non-interacting colloids represent a trivial 
case, it is nevertheless instructive to consider the mean- 
squared-displacement (MSD), characterizing the diffu- 
sive particle motion, both parallel and orthogonal to the 
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flow direction in simple shear [55|- In both the vorticity 
and shear gradient directions, flow has no influence and 
the equihbrium result is recovered, Sz'^ = Sy'^ ~ ^Dq t, 
with Do the single particle diffusion coefficient. In the 
flow direction the MSD is enhanced by a coupling be- 
tween Brownian motion and affine advection, yielding 
6x^ = 2Dot{l + where 7 is the shear rate. The 

physical origin of this enhanced diffusion, termed 'Taylor 
dispersion' [50| , is that the random motion of a given col- 
loid leads to its displacement into planes of laminar flow 
with a velocity different to that of the original point. This 
constant and random 'changing of lanes' leads, on the av- 
erage, to a dramatically increased rate of diffusion in the 
direction of flow. The accelerated rate of mixing achieved 
by stirring a dilute dispersion is thus almost entirely at- 
tributable to local Taylor dispersion. We note also that 
analogous effects arising from flow-diffusion coupling can 
also be identifled in other flow geometries, such as the 
practically relevant case of Poiseuille flow along a cylin- 
drical tube [11] . 



3.2. Dimensionless parameters 

The Smoluchowski equation describes the dynamics 
of spherical colloidal particles dispersed in an incom- 
pressible Newtonian fluid and provides the fundamental 
starting point for all theoretical work to be described 
in the following sections. An appropriate dimensionless 
Reynolds number governing the solvent flow may be de- 
fined as Re = p^R'^/r], with 7 a characteristic fiow rate, 
p is the density, 77 the solvent viscosity and R the col- 
loidal length scale. Due to the small size of the colloidal 
particles Re remains small for all situations of physical 
relevance and the Stokes equations, rather than the more 
complicated Navier-Stokes equations, may thus be em- 
ployed in treating the solvent flow. 

Given that Re remains small, two dimensionless pa- 
rameters are of particular importance in determining the 
equilibrium and nonequilibrium behaviour. The first of 
these is the colloidal volume fraction (p — AnnR'^/S, with 
number density n and particle radius R. The maximum 
volume fraction achievable for monodisperse spheres is 
0.74 corresponding to an optimally packed face-centered- 
cubic crystal structure. For the purposes of the present 
work we will find it convenient to divide the physical 
range of volume fractions into three subregions: (i) low 
packing, (f> < 0.1, (n) intermediate packing, 0.1 < (f) < 
0.494, and (ni) high packing, 0.494 < 0. While this di- 
vision is somewhat arbitrary, it will later prove useful in 
discussing the various theoretical approximation schemes 
currently available. 

The second important dimensionless parameter is the 
Peclet number Pe = 7i?^/2Do The Peclet number 
is a measure of the importance of advection relative to 
Brownian motion and determines the extent to which 
the microstructure is distorted away from equilibrium by 
the flow fleld. In the limit Pe — > Brownian motion 



dominates and the thermodynamic equilibrium state is 
recovered. Conversely, in the strong flow limit, Pe — >■ 
00, solvent mediated hydrodynamic interactions may be 
expected to dominate the particle dynamics, although, in 
practice, surface roughness and other perturbing effects 
turn out to complicate this limit (see section 16.21 for 
more details on this point). 

Finally, we would like to note that there exists a fur- 
ther, nontrivial dimensionless quantity implicit in the 
many-body Smoluchowski equation (|14[) . An increase in 
either the dispersion volume fraction or attractive cou- 
pling between particles is accompanied by an increase 
in the structural relaxation timescale of the system 
characterizing the temporal decay of certain two-point 
autocorrelation functions. This enables the Weissenberg 
number to be deflned as Wi = 7Tq. For intermedi- 
ate and high volume fractions, particularly those close 
to the colloidal glass transition, it is the Weissenberg 
number, rather than the 'bare' Peclet number Pe, which 
dominates certain aspects of the nonlinear rheological re- 
sponse, as has been emphasized in [ssf . For the low vol- 
ume fraction systems to be considered in section [6.21 the 
structural relaxation timescale is set by R^/2Do, leading 
to Pe = Wi. 



3.3. Neglecting solvent hydrodynamics 

In many approximate theories aiming to describe in- 
termediate and high volume-fraction dispersions the in- 
fluence of solvent hydrodynamics beyond trivial advec- 
tion is neglected from the outset. For certain situations 
(e.g. glasses) this approximation is partially motivated 
by physical intuition, however, in most cases, the omis- 
sion of solvent hydrodynamics is an undesirable but un- 
avoidable compromise made in order to achieve tractable 
closed expressions. Accordingly, the expression for the 
probability flux (|15p is approximated in two places, which 
we will now discuss in turn. 

The first approximation is to set Dij — DoSij, thus 
neglecting the influence of the conflguration of the N 
colloidal particles on the mobility of a given particle. For 
low and intermediate volume fraction fluids this may be 
reasonable for Pe <^ 1 but can be expected to break down 
for Pe > 1 as hydrodynamics becomes increasingly im- 
portant in determining the particle trajectories. In par- 
ticular, the near fleld lubrication forces [48| which reduce 
the mobility when the surfaces of two particles approach 
contact play an important role in strong flow and are re- 
sponsible for driving cluster formation and shear thicken- 
ing [5^ (see section [53)) . For dense colloidal suspensions 
close to a glass transition the role of hydrodynamics is 
less clear. For certain situations of interest (e.g. glasses 
close to yield) the relevant value of Pe is very small and 
suggests that hydrodynamic couplings should not be of 
primary importance. 

The second common approximation to (|15p arising 
from the neglect of solvent hydrodynamics is the as- 
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sumption of a translationally invariant linear flow pro- 
file v(r, = K{t)-r, where K{t) is the (traceless) time- 
dependent velocity gradient tensor introduced in Eq.(l8]). 
In an exact calculation the solvent flow field follows from 
solution of Stokes equations with the surfaces of the N 
colloidal particles in a given configuration providing the 
boundary conditions (essentially what is done in Stoke- 
sian dynamics simulation [42,, 51j ) . By replacing this 
solvent velocity field with the affine flow, we neglect the 
need for the solvent to flow around the particles and are 
thus able to fully specify the solvent flow profile from 
the outset, without requiring that this be determined as 
part of a self consistent calculation [g^I • If necessary, the 
assumption of purely affine fiow could be corrected to 
first order. For example, under simple shear fiow the sol- 
vent fiow profile around a single spherical particle is well 
known [31| and could form the basis of a superposition- 
type approximation to the full fiuctuating velocity field. 

It is important to note that the assumption of a trans- 
lationally invariant velocity gradient K{t) is potentially 
rather severe as it excludes from the outset the possibility 
of inhomogeneous flow, as observed in shear banded and 
shear localized states. While physically reasonable for 
low and intermediate density colloidal fluids, the assump- 
tion of homogeneity could become questionable when 
considering the flow response of dynamically arrested 
states, for which brittle fracture may preclude plastic 
flow [6l[. Moreover, it is implicit in the approximation 
\-{r,t) = K{t)-r that the imposed flow proflle acts instan- 
taneously throughout the system. In experiments where 
strain or stress are applied at the sample boundaries a 
flnite time is required for transverse momentum diffusion 
to establish the velocity fleld. Nevertheless, experiments 
and simulations of the transition from equilibrium to ho- 
mogeneous steady state flow have shown that a linear 
velocity profile is established long before the steady state 
regime is approached, thus suggesting that the assump- 
tion of an instantaneous translationally invariant fiow is 
acceptable for certain colloidal systems (65| . 



is clearly inadequate for nonpotential flows (e.g. simple 
shear), for potential flows (e.g. planar elongation) it is 
perhaps tempting to assume such a distribution by cm- 
ploying an effective 'flow potential' Uf (see e.g. [66|). 
The fundamental error of assuming an 'effective equilib- 
rium' description of nonequilibrium states is made very 
clear by the non-normalizability of the assumed distribu- 
tion ^ cxp{—/3{Un + Uf)). These considerations serve 
to emphasize the fact that the only true way to deter- 
mine the distribution function for systems under flow is 
to solve the Smoluchowski equation ([H]). 

For much of the present work we will focus on the re- 
sponse of colloidal dispersions to steady flows. While 
experiment and simulation clearly demonstrate that well 
deflncd steady states may be achieved following a period 
of transient relaxation, it is interesting to note that there 
exists no mathematical proof of a Boltzmann H-thcorcm 
for Eq. (|14[) which would guarantee a unique long-time 
solution for the distribution function. The absence of a 
H-theorem for colloidal dispersions under steady flow is 
a consequence of the hard repulsive core of the particles 
which invalidates the standard methods of proof gener- 
ally applied to Fokkcr-Planck-type equations [13, 113] . 

A further nontrivial aspect of Eq.([T4]) emerges when 
considering the translational invariance properties of the 
time dependent distribution function 5'(t) = 5'(t, {ri}), 
achieved by shifting all particle coordinates by a con- 
stant vector r[ ~ + a (see subsection 7.3 for more 
details). For an arbitrary incompressible flow it has been 
proven that a translationally invariant initial distribu- 
tion function leads to a translationally invariant, but 
anisotropic distribution function ^( t), despite the fact 
that the Smoluchowski operator |15| generating the dy- 
namics is itself not translationally invariant [l8| . Al- 
though the proof outlined in [l^ omitted hydrodynamic 
interactions, it may be expected that the same result 
holds in the presence of hydrodynamics due to the de- 
pendence of the diffusion tensors on relative particle co- 
ordinates. 



3.4. Nonequilibrium states 

In equilibrium, the principle of detailed balance as- 
serts that the microscopic probability flux vanishes, ji = 
^'9i(ln^' -I- PUn) = 0, where Um is the total interparti- 
cle potential energy. This balance between conservative 
and Brownian forces thus yields the familiar Boltmann- 
Gibbs distribution ^I'e = exjp{— PU m) I Z jsi , where Z]^ is 
the conflgurational part of the canonical partition func- 
tion. In the presence of flow (/«(t) ^ 0) there exists a 
finite probability current which breaks the time rever- 
sal symmetry of the equilibrium state and detailed bal- 
ance no longer applies. A nonvanishing probability cur- 
rent thus serves to distinguish between equilibrium and 
nonequilibrium solutions of (|14p and rules out the possi- 
bility of a Boltzmann-Gibbs form for the nonequilibrium 
distribution. While such a Boltzmann-Gibbs distribution 



4. QUIESCENT STATES 

4.1. Hard-spheres 

Theoretical and simulation studies based on Eq. p^ 
have focused largely on the hard-sphere model. In ad- 
dition to being mathematically convenient, the focus on 
this simple model is motivated largely by the availability 
of well characterized hard-sphere-like experimental col- 
loidal systems Q. In the absence of fiow, a system of 
monodispcrse hard-sphere colloids remain in a disordered 
fiuid phase up to a volume fraction of (/) = 0.494, beyond 
which they undergo a first-order phase transition to a 
sohd phase of = 0.545 with face-centered-cubic order 
(see Figure [3]). This unexpected, entropically driven, or- 
dering transition was first observed using molecular dy- 
namics computer simulation in the late 1950's (67| and 
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FIG. 3: A schematic illustration of the phase diagram of 
hard-spheres as a function of volume fraction. Monodisperse 
systems undergo a freezing transition to an FCC crystal with 
coexisting densities (p = 0.494 and 0.545. Polydispersity su- 
presses the freezing transition resulting in a glass transition 
at ~ 0.58, which lies below the random close packing value 
of <;/> ~ 0.64. 



remains a current topic of both experimental and theo- 
retical research (for a recent review see [68|). 

Making the system slightly polydisperse frustrates 
crystalline ordering and suppresses the freezing transi- 
tion. In sufficiently polydisperse systems [69j a disor- 
dered fluid remains the equilibrium state up to a volume 
fraction (j) « 0.58, at which point the dynamics becomes 
arrested and a colloidal glass state is formed. This dy- 
namical transition to a non-ergodic solid is character- 
ized by a non-decaying intermediate scattering function 
at long times for which dynamic light scattering results 
0, @ are well described by the mode coupling theory 
(MCT) [ni- The standard quiescent MCT consists of 
a nonlinear integro-differential equation for the transient 
density correlator which exhibits a bifurcation, identified 
as a dynamic glass transition, for certain values of the 
system parameters [ll| . One of the appealing aspects of 
MCT is the absence of adjustable parameters: All infor- 
mation regarding both the particle interaction potential 
and thermodynamic state point enter via the static struc- 
ture factor, which is assumed to be available from either 
independent measurements or equilibrium statisitical me- 
chanical calculations. For monodisperse hard-spheres, 
MCT predicts a dynamic glass transition at </> ~ 0.516 
when the Percus-Yevick jla | approximation is used to 
generate the structure factor, although other values may 
be obtained using either alternative theories, simulation 
or experiment to determine the static equilibrium struc- 
ture [201 • We note that using MCT together with Percus- 
Yevick structure factors enables a glass transition to be 
studied for monodisperse hard-spheres at volume frac- 
tions above freezing. Neither MCT nor PY theory is 
capable of incorporating crystalline ordering effects and 
both implicitly assume an amorphous microstructure. 

A shortcoming of the quiescent MCT is that it predicts 
an idealized glass transition with a divergent structural 
relaxation time and does not incorporate the activated 
processes which in experiment and simulation studies arc 



found to truncate the divergence. While extensions of 
MCT aiming to incorporate additional relaxation chan- 
nels have been proposed [7l|, [t^] , the underlying micro- 
scopic mechanisms remain unclear. Despite its mean- 
field character, the MCT does capture some aspects of 
the heterogeneous dynamics |73l476| w hich have been ob- 
served using confocal microscopy |77| . 

Finally, we note that a similar scenario of crystalliza- 
tion and dynamical arrest may be observed also in two- 
dimensional systems [1^ [ill . Despite the reduced dimen- 
sionality and new physical mechanisms associated with 
melting in two-dimensions (where the hexatic phase plays 
an important role) the phase diagram for both monodis- 
perse and polydisperse hard-disc systems is qualitatively 
identical to the three-dimensional case illustrated in Fig- 
ure [31 The close analogies between two- and three- 
dimensional systems may be exploited when considering 
nonequilibrium situations for which numerical calcula- 
tions in 3D may prove prohibitively time consuming [79} . 
Viewing a binary mixture as the simplest form of polydis- 
persity, MCT has been employed to study the influence 
of 'mixing' (variations in composition and size ratio) on 
the glass transition of three dimensional hard-sphere [soj 
and two-dimensional hard-disc [8l| systems. These stud- 
ies have revealed intriguing connections between glassy 
arrest and random close packing. 



4.2. Attractive spheres 

The addition of an attractive component to the hard- 
sphere potential can lead to an alternative form of dy- 
namical arrest to either a gel at low volume fraction 
[S^ Issj or an attractive glass state at higher volume 
fractions (s^ . Issj when the interparticle attraction be- 
comes sufficiently strong. The origins of the attractive 
interaction are various e.g. van der Waals forces @ 
or the depletion effect when nonadsorbing polymer is 
added to a dispersion [sol - issj . This form of dynamical 
arrest has been investigated experimentally using both 
dynamic light scattering (see e.g. (ssl - fssj ) and confo- 
cal microscopy (see e.g. |89|). There is now compelling 
evidence both from experiment (90| and simulation [9l| 
that for finite densities gellation occurs via a process 
of arrested phase separation and that only for very di- 
lute, strongly attractive, suspensions does this mecha- 
nism cross over to one of diffusion limited aggregation. 

When applied to attractive colloidal systems the MCT 
predicts a nonequilibrium 'phase diagram' which is in 
good agreement with the results of experiment and qual- 
itatively describes the phase boundary sepcrating fluid 
from arrested states as a function of volume fraction and 
attraction strength [sj, ISSj. Recent studies of systems 
in which the depletion attraction between particles is 
complemented by the addition of a competing long range 
electrostatic repulsion [92j have revealed a rich and unex- 
pected phase behaviour, including stable inhomogencous 
phases [93| and metastablc arrested states [l^l ■ In addi- 
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FIG. 4: The shear viscosity of an aqueous dispersion of col- 
loidal latex as a function of the externally applied shear stress. 
Data for a range of volume fractions are shown, from dilute 
up to a dense colloidal liquid at = 0.50. Shear thinning is 
evident at intermediate stress values as the viscosity of the 
dispersion decreases due to ordering of the particles by the 
flow. At larger applied stresses, for sufHcently high volume 
fraction, the dispersion shear thickens as hydrodynamic lubri- 
cation forces lead to cluster formation and increased disorder. 
(Figure adapted from [o^ ) 



tion, impressive new developments in colloid chemistry 
have enabled the construction of 'colloidal molecules' in 
which the particle surface is decorated with a prescribed 
number of attractive sites, thus rendering the total inter- 
action potential anisotropic [o^. For a review of these 
more recent developments we refer the reader to [96j . 



5. RHEOLOGICAL PHENOMENOLOGY 

As noted in the introduction, dispersions of spherically 
symmetric colloidal particles exhibit a diverse range of 
response to externally applied flow. Much, although not 
all (93, of the generic rheological behaviour of colloidal 
dispersions is captured by the hard-sphere model intro- 
duced in subsection 14. II In order to focus the discussion 
we will consider the special case of hard-spheres subject 
to a steady shear flow. In FigU] we show the results of 
stress controlled experiments performed on a dispersion 
of spherical latex particles dispersed in water at various 
volume fractions, ranging from a dilute 'colloidal gas' up 
to (f) — 0.5, corresponding to a dense colloidal liquid state 
close to the freezing transition [ol] . We note that for the 
experimental steady shear flow data shown in Fig|4] it is 
not significant that the shear stress is employed as the 
control parameter dictating the flow. The quiescent sys- 
tem is ergodic at all considered state points and qualita- 
tively identical results may thus be expected in an analo- 
gous strain controlled experiment, provided that the flow 
remains homogeneous. 



5.1. Zero-shear viscosity 

For each of the volume fractions shown in Fig|4] the 
shear viscosity rj = (Jxy/j is constant for small applied 
stresses (corresponding to small shear rates) and defines 
the zero-shear viscosity 770. The data shown in Fig|4] 
clearly demonstrates that the addition of colloidal parti- 
cles leads to a dramatic increase of 770 above that of the 
pure solvent (note the logarithmic scale in FigU]). 

From a theoretical perspective, there are two alterna- 
tive ways to understand the increase of 770 as a function of 
(j). The first is to relate the viscosity to the flow distorted 
pair correlation functions in the limit of vanishing flow 
rate (sec section [6.ip . The leading order anisotropy of 
g{r, Pe — )• 0) captures the perturbing effect of weak flow 
on the microstructure and thus describes the increase of 
77 in terms of temporally local and physically intuitive 
correlation functions. The second method, referred to as 
either the 'time correlation' or 'Green-Kubo' approach, 
provides an equally rigorous method in which the viscos- 
ity is expressed as a time integral over a transverse stress 
autocorrelation function (see section [7]). Although the 
two approaches are formally equivalent, it is the latter 
which enables a direct connection to be made between rjo 
and the timcscalc describing the collective relaxation of 
the microstructure. 

Within the Green-Kubo formalism the thermodynamic 
coll oida l contribution to the zero-shear viscosity is given 
by [m 
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^= / dtGeq(i), 
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(17) 



where the equilibrium shear modulus is formally defined 
as a stress autocorrelation function 

GcqW = ^^(a,,e"^.*a,,), (18) 

where V is the system volume an d ill ^ is the equilibrium 
adjoint Smoluchowski operator [15|. The fluctuating 
stress tensor element is given by a^y = — J2i ^^''^J ^ a-iid 
the average is taken using the equilibrium Boltzmann- 
Gibbs distribution. 

Eq. lflT)) is an exact Green-Kubo relation which ex- 
presses a linear transport coefficient, in this case the 
shear viscosity, as an integral over a microscopic auto- 
correlation function. For dense colloidal dispersions the 
shear modulus starts from a well defined initial value 
[99! from which it rapidly decays on a time scale set by 
cP/Dq to a plateau. For much later times the modulus 
decays further from the plateau to zero, thus identifying 
the timescale of structural relaxation Tq (see FiglSj. The 
'two step' decay of the time dependent shear modulus is 
a generic feature of interacting systems exhibiting both 
a rapid microscopic dynamics and a slower, interaction 
induced, structural relaxation and is familiar from exper- 
iments and simulations of both colloidal and polymeric 
systems (where the Fourier transform G*{uj) is typically 
considered, rather than G{t) directly). 



11 



Within the idealized mode-coupling theory (MCTVthe 
equilibrium shear modulus (|18p is approximated by jlOlj 



(19) 



where T is the temperature, Sk and S'f. are the static 
structure factor and its derivative, respectively, and 3>fe (t) 
is the transient density correlator defined by 



<i>fc(i) 



Sk 



(20) 



where pk = 6xp(jk • rj). The collective coordinates 
Pk are the central quantity within mode-coupling ap- 
proaches and their autocorrelation (PH)) describes the 
temporal decay of density fluctuations which slow and 
ultimately arrest as the glass transition is approached. 
The mode-coupling approximation (|19p arises from pro- 
jection of the dynamics onto density-pair modes and thus 
expresses the relaxation of stress fiuctuations in terms 
of density fluctuations. Within MCT the correlator is 
approximated by the solution of a nonlinear integro- 
diffcrential equation 



0, (21) 



q'^/Sq and the memory function is 



where Tq 

quadratic functional of ^k(t) which depends upon both 
volume fraction and the static structure factor (which 
serves as proxy for the pair interactions) . Explicit expres- 
sions for the quiescent memory function may be found in 
jl03| . Equation (^1]) predicts that Tq diverges at the glass 
transition volume fraction which then leads, via Eqs. ([T7|) 
and pop , to a corresponding divergence of rjo- While in 
many cases the quiescent MCT give a good account of 
experimental data [12[ the precise nature and location of 
this appa rent divergence remains a matter of debate (see 
e.g. [lOll). 



5.2. Shear thinning 

Turning again to FigU] it is evident that for a given 
volume fraction the viscosity decreases as a function of 
shear rate. This shear thinning behaviour typically sets 
in when the shear rate begins to exceed the inverse of 
the timescale governing structural relaxation, that is for 
values of the Weissenberg number Wi = jTa > 1, where 
7 is the characteristic rate of strain and Tq, is the struc- 
tural relaxation time. Within the range < Wi < 1 the 
system is within the linear response regime and the flow 
rate is sufficently slow that the collective relaxation of 
the microstructure, characterized by, e.g. the decay of 
the transient density correlator ([20]) . is not influenced. 

For Wi > 1 the rate of structural relaxation is en- 
hanced by the flow field. The modulus thus becomes a 
function of the shear rate and the viscosity shear thins. 




FIG. 5: The generalized shear modulus G{t) = G{t, Pe) at 
a vohime fraction relative to the glass transition 4> ~ 4'a ~ 
— 1.16x 10""^ calculated using the extended mode-coupling ap- 
proach [sH]. The final relaxation of the equilibrium modulus 
(blue curve) serves to define the alpha relaxation timescale. 
As the applied shear rate is increased the relaxation timescale 
is reduced. Curves are shown for Pe = 5.5 x lO"'^ (green). 
1.1 X 10"^ (red) and 5.5 x lO"'* (black). The plateau value 
which develops for volume fractions approaching the glass 
transition is indicated by the broken line. The inset shows 
the same data as a function of strain "yt. (Reproduced from 



To incorporate this nonlinear response Eq. (|17p may be 



generalized to 



ri{Pe) 



(Pe) 
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= / dtG{t,Pe), 
Jo 



(22) 



where the functional dependence on Pe has been made 
explicit. The nonlinear modulus is thus defined as 



G{t,Pe) 



k XV ^ "xyh 



(23) 



where $7^ is the adjoint Smoluchowski operator generat- 
ing the particle dynamics [isj . Despite the equilibrium 
averaging employed in (fT8)) . it is important to note that 
an initial stress fiuctuation Cxy evolves to a fluctuation at 
later time t under the full dynamics, including the effects 
of flow. This serves to distinguish the transient stress 
correlator (fT8)) from that which would naturally be mea- 
sured in a computer simulation, where all averaging is 
performed with respect to the full nonequilibrium distri- 
bution function. In the absence of flow 51 = and the 



equilibrium result (|18p is recovered. 

Recent generalizations of the mode-coupling theory 
[2l] | provide approximate expressions for the nonlinear 
modulus G'(t, Pe), see section[71 These approaches incor- 
porate the effects of shear flow into the memory kernel 
responsible for slow structural relaxation and describe 
the speeding up of the rclaxational dynamics. In FigjS] 
we show G{t,Pe) calculated using the generalized MCT 
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[65| for a volume fractions close to (but below) the glass 
transition at various values of the shear rate. For Wi < 1 
the equilibrium result (blue curve) is not influenced by 
the flow. However, for Wi > 1 the longest relaxation 
time becomes dictated by the flow and Tq, ^ 7~^, as is 
demonstrated by the inset to Fig H] which shows the same 
data as a function of strain. From Eg. ([2^ it is clear 
that within this generalized Green-Kubo approach the 
decrease of Tq, with increasing 7 results in shear thinning 
of the viscosity. 

An alternative, although equally valid, viewpoint is 
provided by approaches focusing on the flow distorted 
pair correlations. Exact results for low volume fraction 
dispersions based on the pair Smoluchowski equation (see 
section [6TT|) have shown that shear thinning results from a 
decrease in the Brownian contribution to the shear stress 
[s?! . In dispersions at higher volume fraction the reduc- 
tion in the Brownian stress is manifest in an ordering of 
the particles in the direction of flow which serves to re- 
duce the frequency of particle collisions. Within the shear 
thinning regime evidence for layered or string-like order- 
ing in intermediate volume fraction systems has been 
provided by Brownian dynamics simulations [l04l Il05j | 
and, albeit with different characteristics, by Stokesian 
dynamics s imulation s which include hydrodynamic inter- 
actions [50l . ll06l4l08| . Interestingly, simulations have also 
shown that the flow induced order continues to develop 
following its initial onset. This 'ripening' of the ordered 
phase leads to a time-d ependence of the viscosity known 
as thixotropy and complicates the determina- 

tion of flow curves in both simulation and experiment. 
For each shear rate the measurement time must be suffi- 
ciently long that the viscosity saturates to a plateau value 
before the shear rate is updated. 

We note that the same ordering mechanism discussed 
here for colloidal dispersions would also lead to an analo- 
gous shear thinning scenario for atomic liquids (e.g. liq- 
uid argon). In this case, however, the shear rates required 
to observe such non-Newtonian rheology are several or- 
ders of magnitude larger than those readily accessible in 
experiment. For this reason, non-Newtonian effects in 
atomic systems remain largely a matter of academic in- 
terest. 



5.3. Shear thickening 

Following the regime of shear thinning, a second New- 
tonian plateau develops for which the viscosity attains an 
approximately constant value as a function of shear rate 
and the flow induced ordering of the system continues to 
develop. At higher shear rates the viscosity undergoes a 
rapid increase once a critical value of the shear stress is 
exceeded. Such shear thickening b ehav iour can be either 
continuous [s^ or discontinuous |lll| in character and 
is somewhat counterintuitive in light of the discussion 
presented in 15.21 regarding flow induced microstructural 
ordering and its connection to shear thinning. Suspen- 



sions of nonaggregating particles at intermediate volume 
fractions generally show reversible shear thickening, how- 
ever the details of the increase in viscosity depend upon 
the details of the system (particle type, solvent etc.) as 
well as the thermodynamic control parameters |ll2| . 

In section [5?^ we noted that the onset of shear thinning 
occurs for values of the Weissenberg number Wi > 1, 
reflecting the essential competition between flow and 
structural relaxation. In contrast, the onset of shear- 
thickening behaviour is determined by the value of the 
bare Peclet number Pe, thus serving to highlight the dif- 
ferent mechanisms dominating the physics of thinning 
and thickening states. In experiment, the difference in 
scaling of Wi and Pe with particle diameter d enables the 
extent of the second Newtonian plateau separating shear 
thinning and shear thickening regimes to be controlled as 
a function of particle size. When the microstructure un- 
der flow is also of interest, the contraints of instrumental 
resolution (in e.g. confocal microscopy) place additional 
limits on the particle size which have also to be taken 
into consideration. 

The earliest theoretical explanations of shear thicken- 
ing in colloidal dispersions proposed that the observed 
viscosity increase i s the co nsequence of an order-to- 
disorder transition [TT^ . Within this picture, the 
ordered planes of particles which form within the shear 
thinning regime, and which persist throughout the sec- 
ond Newtonian viscosity plateau, begin to interact via 
hydrodynamic coupling at sufficiently high shear rates. 
This interaction pulls particles out of the layers, lead- 
ing to increased particle collisions, disorder, and a conse- 
quent increase in viscosity. The onset of shear thickening 
is thus identifled with the h ydro dynamic instability of a 
layered microstructure (see |ll5j | for a discussion of this 
issue). 

Despite the intuitive appeal of interpreting shear thick- 
ening as an order-to-disorde r transiti on, questions were 
raised by the experiments of |ll5l Ill6l | in which a specific 
system (electrostatically stabilized Latex particles in gly- 
cols) was found to display shear thickening in the absence 
of an ordered phase. These results suggested that while 
an ordered phase may well precede the shear thickening 
regime as the shear rate is increased, it is not a necessary 
prerequisite. According to these findings, shear thicken- 
ing occurs via an independent physical mechanism and 
is not simply related to a loss of microstuctural order. 
Further insight into the microscopic mechanism underly- 
ing shear thi cken i ng w as provided by Stokesian dynamics 
simulations |104 Ill7j | which identified the formation of 
hydrodynamically bound particle clusters at high shear 
rates. Such 'hydroclustcrs' form when the shear flow is 
sufficiently strong that the particle surfaces are driven 
closely together. At such small separations the hydrody- 
namic lubrication forces dramatically reduce the relative 
mobility of the particles such that they remain trapped 
together in a bound orbit (a point which we will later re- 
visit in section 15^ . Transient shear-driven hydroclustcrs 
would appear to be the defining feature of shear thick- 



13 



^o°°o°o^o 



ooooo^o 

ooo Oo 
oqOooo o 

OOOOqO o 



oo^ o oo 



ocP 
o 

o 
o 



o 



oo 



Newtonian 



Shear tliinning Shear thickening 



Pe 



FIG. 6: Schematic iUustration of the microstructural order 
and disorder induced in a dense colloidal dispersion by shear 
flow. At low values of Pe (leftmost configuration) the viscos- 
ity remains constant as diffusion is able to restore the equi- 
librium microstructure more rapidly than the shear flow can 
disrupt it. At intermediate shear rates (central configuration) 
the rate of shear exceeds the rate of structural relaxation, 
Wi > 1 leading to microstructural ordering and shear thin- 
ning. At high shear rates (right configuration) hydrodynamic 
lubrication forces lead to particle clustering which strongly 
enhances the hydrodynamic contribution to the viscosity and 
result in shear thickening. 



ened states and experi ment al evidence for their impor- 
tance is accumulating jll8l |. Nevertheless, opinion re- 
mains divided regarding the fundamental mechanisms at 

work fngt - 

The hydrodynamic mechanisms described above give 
rise to a continuous, albeit rapid, rise in the viscosity as 
a function of shear rate. An alternative scenario may 
arise when, at some critical value of the shear rate, the 
viscosity exhibits a discontinuous jump as the system be- 
comes jammed [H, [HI [TMIlll . While it is anticipated 
that hydrodynamics will be relevant for the description 
of shear thickening at intermediate volume fractions (e.g. 
< (p < 0.5, as considered in FigH} alternative mecha- 
nisms may bec ome imp ortant upon approaching the glass 
transition. In 123l4l25 a 'schematic' mode coupling the- 
ory similar to those to be discussed in section [7] was de- 
veloped, in which a coupling to stress was introduced 
into the nonlinear equations determining the decay of the 
transient density correlator. Upon varying the model pa- 
rameters a range of rheological behaviour was revealed, 
including both continuous and discontinuous shear thick- 
ening, as well as a jamming transition to a non-ergodic 
solid state. Within this picture, shear thickening and 
jamming are viewed as a type of stress induced glass 
transition, for which the applied stress inhibits particle 
motion, even in the absence of hydrodynamics. A num- 
ber of works have suggested a relat i onsh ip between shear 
thickening and jamming 0, 0, llOSl Il26l | although details 
of the connection between hydrodynamic cluster forma- 
tion and jamming transitions o f the kind more familiar 
from studies of granular media |127| remains unclear. 

The addition of an attractive component to the 
strongly repulsive colloidal core can lead to gel formation 
and irreversible flocculation (see section 14. 2p . For such 
systems shear thickening is generally not observed, as the 
increase in the hydrodynamic contribution to the viscos- 
ity with increasing shear rate is more than compensated 
for by the decrease in the thermodynamic contribution 



arising from the attraction (see e.g. [I28l |). The generic 
behaviour of gel and floe states is thu s mo notonic shear 
thinning as a function of shear rate (ll2| . It is there- 
fore surprising that recent experiments using attractive 
carbon black particles jl29l | have identified a rich shear 
thickening behaviour for which the viscosity increases be- 
yond a critical value of the applied shear stress. In this 
case an additional physical mechanism has been proposed 
by which the forces exerted by shear flow cause floes to 
break apart, leading to an increase d sur face area and thus 
greater hydrodynamic dissipation [l29l |. 

As pointed out in section 15.21 shear-thinning is not 
unique to colloidal system and can also be observed, al- 
beit at high shear rates, in simple atomic liquids. In 
contrast, shear thickening of the type discussed above is 
not found in atomic systems (for which the 'solvent' is a 
vacuum) and demonstrates clearly the breakdown of the 
correspondence between colloidal and simple liquids for 
strongly noncquilibrium states. Indeed it is quite clear 
that the view of colloids as 'big atoms' jl3C| will only 
hold in situations for which the influence of the solvent is 
negligable and that new physics may emerge when the 
role of hydrodynamic interactions becomes significant. 
We note that an alternative type of shear thickening has 
been observed at high shear r ates in mo lecular dynamics 
simulations of simple liquids [I3ll4l34 . In these simu- 
lations a profile unbiased thermostat was employed to 
remove artifacts which may arise when a linear flow pro- 
file is assumed. Shear thickening was observed in simu- 
lations performed at constant volume, but not in those 
performed at constant pressure. 

Finally, we would like to note that the onset of shear 
thickening at high flow rates has been associated with 
unexpected behaviour of the first normal stress differ- 
ence A^i = a^x — o-yy ■ Typically, dispersions at low or 
moderate shear rate exhibit a positive value of Ni, indi- 
cating that a Weissenberg (or 'rod climbing') effect would 
be observed in shear experiments performed in a Couette 
geometry (H . Experiments on dense colloi dal dispersions 
with repulsive interactions |l2ll Il35l Il36| have revealed 
that A^^i can change sign from positive to negative upon 
increasing the flow rate into the regime where the viscos- 
ity shear thickens. Similar behavi our h as been observed 
in Stokesian dynamics simulations |l37l | and in numerical 
solut ions of the Smoluchowski equation for dilute systems 
jl38l | . In contrast to these results for purely repulsive in- 
teractions, recent experiments on attractive flocculated 
colloidal dispersions display a monotonicall y in creasing 
A'^i throughout the shear thickening regime |129| . 



5.4. Yield stress 

For colloidal fluid states {0 < (f> < 0.494) the data pre- 
sented in Fig|4] represent the generic phenomenology of 
dispersions of strongly repulsive colloids under shear flow. 
In fact, this behaviour is not limited to simple shear. 
Qualitatively identical behaviour is found in Stokesian 
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FIG. 7: Many commonly studied colloidal particles (e.g. 
PNIPAM) consist of a polymeric core grafted with a layer of 
polymer (schematically represented on the left) which serves 
to stabilize against flocculation. Given a sufficiently dense and 
crosslinked polymer brush the particles exhibit a strongly re- 
pulsive effective potential interaction approximating that of 
hard-spheres with radius R. The ability of the solvent to pene- 
trate into the brush results in a hydrodynamic radius Rh < R- 
These effects may be mimiced by a simple hard-sphere model 
(sketch on the right) in which R/Rh can be used to control 
the influence of hydrodynamic interactions. 



dynamics simulations for the extensional viscosity [23| of 
dispersions of hard-sphere colloids under steady exten- 
sional flow, with the shear rate replaced by the rate of 
Hcncky strain [ssf . 

As already noted in section fS-li increasing the volume 
fraction of a colloidal liquid leads to a strong increase 
in the zero-shear viscosity. Assuming that crystalliza- 
tion has been suppressed by polydispersity, the volume 
fraction can then be further increased, eventually result- 
ing in an apparent divergence of the zero-shear viscosity, 
either at the glass transition vol ume fraction (accord- 
ing to mode-coupling theory [TTI. Il0l| ). or some h ighe r 
volume fraction approaching random close packing |102| . 
The variation of the viscosity as a function of shear rate 
for volume fractions ranging from 0.45 to 0.57 is demon- 
strated in more detail by the data shown in FigIS] These 
experiments were performed on a system of poly (ethylene 
glycol) -g rafte d polystyrene colloidal particles dispersed 
in water |l39l |. For the two lowest volume fractions con- 
sidered ((/) = 0.45 and 0.48) a clear zero-shear viscosity 
may be be identified from the low shear rate plateau, with 
shear thinning evident at higher shear rates for = 0.48. 
As the volume fraction is increased above 0.48 the low 
shear rate plateau moves to smaller rates, out of the 
experimental window of resolution, and the dispersion 
shows shear thinning over the entire range. Analysis of 
the intensity correlation function (related to the tran- 
sient density correlator (j20p ) measured using dynamic 
light scattering leads to an estimate of the glass transi- 
tion for this system of 0.53 < cf>g < 0.55, somewhat lower 
than the typical value (j>g ~ 0.58 obtained for PMMA 
hard-sphere-like colloids. The viscosity data for the two 
highest volume fractions (0 = 0.55 and 0.57) are consis- 
tent with a divergence in the zero-shear viscosity at the 
glass transition, as predicted by the MCT (see section 



According to the extended MCT [HI, [H, [m^, for 

glassy states the slowest relaxation time is Tq, ^ 7""'^, 
which leads, via Eq.dig), to ?7 G{t oo)7"\ where 
G{t — >• oo) is the plateau modulus (see FiglS]), thus re- 
producing the power law decay of the viscosity demon- 
strated by the data in FiglHl For the idealized glassy 
states considered by MCT, where Tq is infinite in the ab- 
sence of flow, this power law dependence extends to the 
limit 7—^0, resulting in a true divergence. In real col- 
loidal experiments, higher order relaxation processes will 
always endow the quiescent system with a finite value of 
Ta and the viscosity divergence will be truncated. The 
low shear divergence of the viscosity and power law shear 
thinning 77 ^ 7"^ suggested by Fig [5] are supported by 
independent experiments performed on thermosensitive 
core-shell particles [63ll63|. However, some recent exper- 
iments on sterically stabilized PMMA particles provide 
contradictory evidence and have suggested a nontrivial 
dependence of the relaxation time on shear r ate, namely 
Ta ^ 7 "'^i which remains to be understood [14ll |. 

The inset to Fig|S] shows the shear stress as a function 
of shear rate and provides an alternative representation of 
the viscosity data shown in the main panel. For the two 
highest volume fraction samples (cj) = 0.55 and 0.57) the 
shear stress becomes constant for the lowest shear rates 
considered, thus identifying a dynamic yield stress for 
glassy states. Both the shear thinning as a function of Pe 
and appearance of a dynamic yield stress as a function of 
(j> evidenced by Fig[5]are well described by the extended 
MCT 

The relationship between the dynamic yield stress and 
the more familiar static yield stress mirrors that be- 
tween stick and slip friction in engineering applications 
l^^stat ^dyn ^j^^g cxpccted). Indeed, it may be 

argued that the dynamic yield stress is, in fact, a more 
well defined quantity than the static yield stress. The lat- 
ter is typically defined as the step stress amplitude which 
must be exceeded such that the system will fiow at long 
times [l42i] . The point of static yield may therefore be de- 
pendent upon details of the system preparation, with the 
consequence that nonstationary properties, such as sam- 
ple age in colloidal gla sses, could influence the outcome 
of a given experiment |l43l |. Moreover, the existence of 
creep motion, for which the strain increases sublinearly 
with time, makes difficult an unambiguous identification 
of the static yield stress. In contrast, the dynamic yield 
stress is defined as the limiting stress within a sequence of 
ergodic, fluidized steady states and is thus independent 
of prior sample history. 

It is apparent from Ea. (|22[) that a dynamic yield 
stress can only exist in the event that ^ 7"", with 
1/ = 1. Values of 1^ less than unity result in a shear stress 
o'xy(7 — >■ 0) = 0, despite the fact that the viscosity di- 
verges. We thus note that a low shear rate divergence of 
the viscosity is a necessary but not sufficient condition for 
th e exis tence of a yield stress. While the results presented 
in [l4l| apparently cast doubts on the existence of a dy- 
namic yield stress for certain colloidal glasses, complica- 
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FIG. 8: Shear thinning and the dynamic yield stress of a con- 
centrated aqueous dispersion of poly(ethylene glycol) -grafted 
polystyrene colloidal particles which, to a good approxima- 
tion, behave as hard spheres. The main figure and inset show 
the viscosity and shear stress, respectively, as a fu nction of 
shear rate. (Reprinted with permission from [l39l |. Copy- 
right: American Physical Society) 



tions due to inhomogeneous, sh ear l ocalized flow make 
this a subject of ongoing debate |144| |. 

The interplay between static and dynamic yield has 
been investigated in simulation studies of a glass form- 
ing binary Lennard- Jones mixture (the Kob - Anderson 
model) using molecular dynamics simulations jl45l Il46l | . 
In these simulations the mixture was confined between 
two atomistic walls, one of which was then subjected to 
either a constant stress or constant strain in order to 
induce shear flow. It is important to note that due to 
the application of shear through the boundaries, the flow 
profile within the confined fiuid/glass is an output of the 
numerical calculation and is not constrained to be linear. 

In Figiniwe show the simulated flow curve for this sys- 
tem (analogous to that shown in the inset to FiglS]) for a 
glassy statepoint, calculated by ap plying a f ixed rate of 
strain to one of the bounding walls jl45l Il46j . When the 
shear stress is plotted as a function of the total strain rate 
(which may differ from the local rate of strain) a dynamic 
yield stress can clearly be identified. It was observed in 
the simulations that at sufficiently low (total) shear rates, 
an inhomogeneous flow profile develops in which a static 
layer coexists with a fiuidized region exhibiting a linear 
flow profile. In a complementary set of simulations a 
lower bound for the static yield stress was identified by 
slowly (stepwise) increasing the shear stress until viscous 
flow could be detected at long times. The static yield 
stress thus obtained was found to provide a criterion for 
determining the onset of inhomogeneous flow. These ob- 
servations may be consistent with experi ment s on PMMA 
colloids exhibiting inhomogeneous flow |l4l| but are ap- 



parently at odds with experiments on core-shell particles 
which do not give indications of banding or shear local- 
ization effects [il, 111, Il39j . While these discrepancies 
remain to be understood, it seems possible that the soft- 
ness of the potenti a.1 int eraction in the core-shell systems 
studied in [63l [6^ Il39l | may play a role in maintaining 
homogeneous flow. 



6. THEORETICAL APPROACHES TO FLUID 
STATES 

There currently exist several alternative theoretical ap- 
proaches to first-principles calculation of the microstruc- 
turc and macroscopic rheology of colloidal dispersions 
subject to externally applied flow. Each of the available 
approximation schemes is tailored to capture the phys- 
ically relevant aspects of the correlated particle motion 
within a restricted range of volume fractions. Theories 
aiming to treat low and intermediate volume fraction dis- 
persions take as their common starting point the pair 
Smoluchowski equation, which is an exact coarse grained 
reduction of the many-body Smoluchowski equation (|14p . 
At high volume fractions close to the glass transition the 
pair Smoluchowski equation no longer provides a conve- 
nient starting point and an alternative approach capable 
of capturing slow structural relaxation is required. This 
is provided by the recently develo ped int egration through 
transients mode-coupling theory |l7l . [isl [2ll . Il40j | . In the 
following, we will first introduce the pair Smoluchowski 
equation before proceeding to follow the 'volume fraction 
axis' to give an overview of the current state of research 
on the theory of fiowing states. 



6.1. Pair Smoluchowski equation 

While Eq. p^ provides a well defined microscopic dy- 
namics, it has been found useful to start from an equiva- 
lent coarse grained level of description by integrating out 
uneccessary degrees of freedom from the outset. Assum- 
ing spatial translational invariance, integration of Eq. (jl4p 
over the center-of-mass coordinate of a pair of particles 
and the remaining N~2 particles leads to an equation for 
the flow distorted pair correlation function as a function 
of r = r2 - ri (see e.g. [57,, ,147nl49il ) 



dgiy) 
dt 



Vr- [v(r)g(r)-D(r)-V,g(r)] 
= -V,- [D(r)./3F(r)g(r)] 



(24) 



where we have suppressed explicit time-dependence in 
the function arguments for notational convenience and 
where we have introduced the gradient operator = 
V2 — Vi- The conditional probability to flnd particles at 
coordinates r3 • • • r , given that the first two are known 
to be at locations ri and r2, respectively, is given by 
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P(r3, • • • ,rAr|ri,r2) = P(ri, • • • , r7v)/P(ri, and is re- 
quired to calculate the functions F(r),D(r) and v(r) en- 
tering Eq.([24|). The first of these functions, F(r), de- 
scribes the force acting between our chosen pair of parti- 
cles due to both direct potential interaction v{r), taken 
here to be pairwise additive, and indirect interactions 
transmitted via the surrounding N~2 particles 

F r = -\/r v{r) - - / rfra ^ — x 

2 J .g(ri,r2) 

X (V2i'(|r2 - ral) - Vi«(|ri - r^D), (25) 

where 5^^^(ri,r2,r3) is the nonequilibrium triplet distri- 
bution function. The diffusion tensor is similarly ob- 
tained by conditional averaging and contains details of 
the hydrodynamic interactions 

D(r)=2i?o(5G(r) + (5-5)iJ(r)), (26) 

where rr denotes a dyadic product and the scalar hydro- 
dynamic functions G(r) and H{r) remain to be specified. 
Finally, the relative velocity of a pair of particles is given 

by 

v(r) =K-r + C(r):K (27) 

where ii is the symmetric ratc-of-strain tensor k = (k + 
k'^)/2 and C is the (third rank) hydrodynamic resistance 
tensor describing the disturbance of the affine flow due 
to the presence of the particles 

C(r):K =-r —Air) + (S--) S(r) 

\ r \ r / r j 

(28) 

The tensor C arises from purely geometrical considera- 
tions and is not material specific. It is interesting to note 
that the addition of hydrodynamic interactions prevents 
advection leading to unphysical hard core overlap. A pair 
of approaching particles thus 'flow around' each other in 
the solvent flow, an effect taken care of by the second 
term in Eq. ((?7)) . 

In the dilute limit, much is known about the hydro- 
dynamic functions A, B, G and iJ, as only an isolated 
pair of spheres must be considered. For both large and 
small separations analytical expressions for these func- 
tions exist [i^ and are supplement ed by tabulated nu- 
merical data for intermediate ranges [l50l | . At higher vol- 
ume fractions approximations are required to obtain the 
hydrodynamic functions and a number of schemes have 
been developed which aim to incor porate the effec ts of 
many-body hydrodynamics H ^ [lo3, [HI [HI. In 
the absence of hydrodynamic interactions A = B = Q 
and G = H = 1 leading to considerable simplification. It 
should be noted that C = in this limit, with the con- 
sequence that affine motion alone can lead to hard-core 
overlaps. While an exact treatment of the thermody- 
namic part of the problem would lend such unphysical 




FIG. 9: Molecular dynamics simulation results for a glass 
forming binary Lennard- Jones mixture. The full line shows 
the shear stress as a function of the total shear rate under 
conditions for which a constant rate of strain is applied to 
one of the bounding walls. The square (labelled ay) indi- 
cates the static yield stress obtained by sequentially increasing 
the applied stress until the system begins to exhibit viscous 
flow. For values of the stress between dynamic and static yield 
points, the system was found to e xhib it inhomogeneous flow. 
(Reprinted with permission from [l45l ]. Copyright: American 
Institute of Physics.) 



configurations zero statistical weight, care must be exer- 
cised in approximate treatments which may satisfy only 
partially this important geometrical constraint. 

Although the coarse grained pair Smoluchowski equa- 
tion p4|) is still exact (under the assumption of homo- 
geneity), it does not provide a closed expression for the 
microstructure, as encoded in ^(r). Evaluation of the 
integral term required to determine the force ([25l) de- 
mands knowledge of the nonequilibrium triplet distribu- 
tion function, which remains unknown and contains the 
residual influence of the surrounding particles which have 
been integrated out. This situation is familiar from the 
BBGKY heirarchy [l^ for which the triplet correlations 
must be approximated in terms of the pair correlations 
(using e.g. the Kirkwood superposition approximation) 
in order to arrive at a closed equation. In recent years, 
accurate approximations for the equilibrium triplet cor- 
relations of certain model systems have been developed 
[l53l Il54l |. Less is known regarding the nature of the 
triplet correl atio ns in nonequilibrium situations. Recent 
simulations |l56l | using accelerated Stokcsian dynamics 
[13, [HI have revealed the existence of aligned particle 
triplets under shear and it may be hoped that such mi- 
crostructural insights will eventially lead to improved 
theories by guiding the development of approximate clo- 
sures for the triplet correlations. Some of the approaches 
to be reported in section 16.31 have attempted to make 
progress in this direction by approximating explicitly the 
integral term on the right hand side of (P^ . 
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We note that, although the assumption of spatial ho- 
mogeneity underlying (j24p is mathematically convenient 
(for a translationally invariant system the only physically 
relevant coordinate is the separation vector r = r2 — ri) 
it may not be appropriate under all conditions. The pres- 
ence of spatial inhomogcncity induced by either external 
potential fields, shear banded or shear localized states 
complicates the coarse graining procedure, resulting in 
an inhomogeneous version of p4)) . While these issues 
should pose no difficulty at low volume fractions, for 
which the right hand side of ([M)) may be disregarded, 
caution should be exercised when treating systems at 
higher volume fraction. 

Once g{r) is known, calculation of the stress tensor 
describing the macroscopic rheological response becomes 
possible. Although exact expressions relating the stress 
tensor to to the flow distorted microstructurc arc known 
formally, the situtation is complicated by the appear- 
ance of unknown conditionally averaged hydodynamic 
functions in the expressions. However, reliable approxi- 
mations for these functions are available and enable the 
stress to be evaluated directly from g{r) . For the sim- 
pler case of a system interacting via a pair potential and 
in the absence of hydrodynamic interactions, the stress 
tensor may be completely determined by a simple inte- 
gral over the pair correlation function |155l | 

Tl'^ f TY 

a = -nkeTl + Sj^^^k dr —v'(r)g(r), (29) 

2 J r 

where v'{r) is the derivative of the pair potential, 77^ the 
solvent viscosity, 1 the identity matrix and rr denotes a 
dyadic product. Note that the second term on the right 
hand side of Eq. assumes that the particles possess a 
well deflncd hard-core from which the solvent is excluded. 

In the absence of flow the stress tensor is diagonal with 
the osmotic pressure given by 11 = — Trcr/3. For a sys- 
tem of pure hard-spheres the familiar equation of state 
pn/n = 1 + A(j)g{d) is thus recovered. When under shear 
flow Ea. (pn| yields a shear viscosity due to the colloids 
77 = 5 (f)ris/2 + 0(41'^), where the first term corresponds 
to Einstcins classic dilute limit result |157l | and the cor- 
rections to higher order in (f) come from the anisotropy of 
g{r) inside the integral term. It is clear from Eg. ((29)) that 
flow induced microstructural anisotropy can give rise to 
the finite normal stress differences Ni = a^x ~ Cyy and 
— Cyy" CTzz characteristic of non-Newtonian rheology. 
The dyadic weight factor entering the integral term has 
the consequence that if g{r) possesses a mirror symmetry 
about the 2; = plane then the integral term will be equal 
to zero and the rheology will thus be Newtonian. While 
this 'forc-aft' symmetry of the pair distribution function 
is an exact mathematical consequence of the 'pure hy- 
drodynamic limit', in which the moti on of the particles 
is determined by Stokes flow alone jl58l Il59l |. chaotic 
many-body particle motion and experimental perturba- 
tions, such as particle surface roughness, present in real 
colloidal systems break the symmetry and result in a non- 
Newtonian rheology [s^ . 



6.2. Low volume fraction 

Efforts to obtain a microscopic understanding of colloid 
rheology began with the seminal 1906 work of Einstein in 
which it was shown how the shear viscosity of a dilute dis- 
persion of hard spherical colloids increases with colloidal 
volume fraction, assuming that both the volume frac- 
tion and the shear rate remain small {rj = 775(1 -I- 5(/)/2)) 



157| . Einstein's study addressed the one-body problem 



of a single colloid suspended in a Newtonian fluid. The 
next step is naturally to consider the interaction between 
pairs of colloidal particles, thus making possible a discus- 
sion of the pair correlation functions and their relation 
to rheological functions at low volume fraction. Study of 
the two-particle dilute limit was initiated by Batchelor 
P [l5i4l60t whose fundamental work formed the basis 
for the more recent investigations by Brady and cowork- 
ers (57L [l38l . [l6H . 

At low volume fraction, Eg. ip^ admits anal ytica l so- 
lution in the limits Pe and Pe 00 [stI . I161| | and 
precise numerical results have been obtained for interme- 
diate values of Pe, bot h with and without hydrodynamic 
interactions jl38l Il64l | . Exact results are made possible 
by the fact that, in the dilute limit, triplet correlations in 
the pair Smoluchowski equation may be neglected leading 
to a closed expression for g{r). Neglecting the difficult 
integral term in Eg.ip^ incurs an 0{(j)) error which be- 
comes irrelevant as </>—>■ and yields a closed equation for 
g{r) which is exact to lowest order in the volume fraction 
and valid for all Pe values. For the simple special case 
of hard-spheres Eg. ([M)) thus reduces to the equation of 
motion 



dgjr) 
dt 



Vr ■ [v(r)g(r) - Pe-^D{r) ■ V.g(r)] = 0, (30) 



where wc have scaled distance and time with particle ra- 
dius and flow rate, respectively, such that Pe appears 
explicitly. In order to fully specify the problem Eg. ([50)) 
must be supplemented with appropriate boundary con- 
ditions enforcing both the requirement that the particles 
do not penetrate, via a no-flux condition at r = d, and 
that g{r) ^ 1 as r 00 Q. The first of these bound- 
ary conditions is clearly an exact physical requirement 
and is valid also at higher volume fractions. The second 
condition assumes the decay of 'wake' structures which 
develop in 17 (r) downstream from the reference particle 
at higher flow rates. Detailed analysis of ([30]) has shown 
that the range of the wake scales linearly with Pe, thus 
justifying the choice of boundary conditions. 

For systems interacting via a spherically symmetric 
pair potential it can be shown that, regardless of vol- 
ume fraction, in the weak flow limit Pe — >■ a steady flow 
field acts on the spherically symmetric equilibrium distri- 
bution gcq{r) to produce an 0{Pe) perturbation j6l. ll60l | 



.g(r) = gcq{r) 



1 



(31) 



where k = n/y/ 2k : k, with k, is defined below Eq. 
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FIG. 10: For mathematically perfect hard-spheres with a 
hydrodynamic radius equal to the radius of potential interac- 
tion, Rh = R, particle pairs exhibit well defined trajectories. 
Taking a reference frame in which one particle is fixed at the 
origin (yellow) the second particle (red) follows the trajecto- 
ries shown in the 'pure hydrodynamic limit' of large Pe [Slj . 
The figure shows some sample trajectories in the z = flow- 
gradient plane (see e.g. [^). Closed orbits are indicated in 
red and open in black. The apparent fore-aft mirror symme- 
try gives rise to Newtonian rheology. 



In the dilute limit 5'oq('') = 0{—r — 1) and substitution of 



nj into ([5n| yields a differential equation for the dimen- 
sionless funetion f{r) which has been solved for several 
interaction potentials of interest Q . Although analytical 
expressions exist for certain systems, a numerical inte- 
gration is still required to obtain /(r) in the specia l cas e 
of hard-spheres with hydrodynamic interactions ll6Cll | . 

The extension of Eq. ([5T|) to higher order in the Peclet 
number has been analyzed in considerable detail |l6l| . 
It can be shown that g{r) has a regular perturbation 
expansion to O(Pe^) but that calculation of the next or- 
der term requires singular perturbation theory, yielding 
an 0{Pe^^'^) correction. The calculation of higher or- 
der terms in the Pe expansion requires use of matched 
asymptotic expansions which rapidly become intractable 
and make preferable a numerical solution of (13pp. The 
expansion of the distorted structure is given by (16ll | 



g{r) = 1 + hPe + hPe^ + h/2Pe^'^ 



(32) 



where comparison with (j3ip . and noting that gcqir > 
2R) = 1 for hard-spheres at low volume fraction, enables 
identification of the coefficient /i. To 0{Pe) the rhe- 
ology is predicted to be Newtonian with normal stress 
differences identically equal to zero [3l| . Non- Newtonian 
rheology first occurs at 0{Pe^), which is sufficient to 
capture both non-zero normal stresses and the first flow 
induced correction to the osmotic pressure |l6l| . 

Analytic solutions to Eg. ([50)) exist also in the 'pure 
hydrodynamic limit' of strong flows (Pe^oo) and have 
highlighted the subtle balance between hydrodynamic 
and poten tial forces in determining the rheological re- 
sponse [sll, Il59| . For large Peclet number steady flows the 
solution of Eq. (|3D]) is well approximated by the solution 



of the the simplified equation V • [v(r) g{r)] = (subject 
to the boundary condition g(r) = 1 at r = oo), despite 
the fact that the approximation neglects the boundary 
layer and thus violates the no-flux condition at contact. 
Subject to certain conditions on the trajectories of par- 
ticle pairs, Batchelor and Green proved the surprising 
result that the simplified equation predicts a spherically 
symmetric radial distribution f uncti on for hard spheres, 
leading to Newtonian rheology |l59l |. This clear predic- 
tion is a direct consequence of the fore-aft symmetry of 
.g(r) (see subsection 16. ip inherent in the assumed Stoke- 
sian solvent flow. In Figure lTUl we show sample trajec- 
tories of a hard-sphere in shear flow as it moves around 
a second sphere held fixed at the coordinate origin. Of 
particular interest is the existence of closed trajectories 
along which the particles become trapped in bound orbits 
and which are connected to the lubrication force required 
to displace solvent from the region between the particles 
[48| . In fact, the lubrication force acting between a pair 
of perfect spheres at separation r shows a divergence. 
Flub ~ ('^/-R ~ corresponding to surface contact. 

Crucially, the time-reversibility of the Stokes equations 
dictating the solvent flow implies that the force required 
to push particles togther is identical to that required to 
pull them apart, with the consequence that particle tra- 
jectories exhibit the fore-aft mirror symmetry apparent 
in FigHni 

De spite the sound mathematical evidence provided in 
jl59j |. serious doubts were cast by subsequent experi- 
ments on intermediate volume fraction hard-sphere-like 
colloidal dispersions, which seemed to contradict the the- 
oretical predictions by id entify ing a non-Newtonian rhe- 
ology at large flow rates [l62| . It should be noted that 
the reversibility of Stokes flow implies that fore-aft sym- 
metry in the pure hydrodynamic limit holds also for finite 
volume fractions an d so the value (j) = 0-4 employed in 
the experiments of jl62| cannot be held responsible for 
the apparent discrepancy. The situation was eventually 
resolved by Brady and Morris, who analytically identified 
a boundary-layer in the region close to particle contact 
in which Brownian motion balances advection 1571. The 



analysis of [57[ indeed recovers the findings of |159l | m 
the case that the hydrodynamic radius is equal to the 
excluded volume radius, as would be the case for mathe- 
matically perfect hard-spheres with no surface roughness. 
However, when the excluded volume radius exceeds the 
hydrodynamic radius, even by a very small amount, the 
residual Brownian motion within the anisotropic bound- 
ary layer of g{r) leads to a non-Newtonian rheology in 
the strong flow limit. FiglT] shows a sketch of the model 
employed in (STj . The physical origin of these symmetry 
breaking surface effects remains an open problem and is 
likely to be a function of various system specific parame- 
ters (e.g. surface roughness). We note that the existence 
of a boundary-layer structure was originally identified in 
studie s of the distorted structure factor of colloids under 
shear Il63ll. 



The analytical results for Pe — > and Pe— > oo for hard- 
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FIG. 11: The pair contributions to the relative shear vis- 
cosity of a dilute colloidal dispersion under steady shear flow 
for three values of the ratio R/Rh (see Fig[7]). To second 
order in the relative viscosity is given by rjr = rj/rjs = 
l + 5<j)H/2 + A{Pe,R/RH)<p'^, where 0h = {RH/Rf(t>, and 
A{Pe, R/Rh) can be divided into hydrodynamic and Brown- 
ian contributions A — Ah + Ab- As R/Rh 1 the dispersion 
shows shear thickening at large Pe values due to the increase 
in the hydrodynamic contribution to the shear stress. As the 
value of R/Rh is increased the increase in Ah becomes bal- 
anced by the decrease in and only shear thinning remains. 
(Figure adapted from [l38l |l 



spheres under steady shear have been both confirmed and 
sup plemented by full numerical solutions at all values of 
These accurate numerical studies revealed 
that the dilute dispersions described by Eq. (pO| demon- 
strate not only shear thinning and finite normal stresses 
at intermediate flow rates (accessible from the expansion 
(|32|) t o O jPe^)), but also shear thickening at high flow 
rates [l38l | . Shear thinning in dilute systems is correlated 
with a nonvanishing distortion of the structu re fac tor in 
the plane perpendicular to the flow direction jl65| . 

In Figure [11] we s how some of the numerical results 
obtained in Ref. |l38l | for the shear viscosity of a dilute 
dispersion as a function of Pe. In these calculations the 
effective sphere model sketched in FigIT] was employed 
and results are shown for three values of the ratio of po- 
tential to hydrodynamic radius, R/Rh- To second order 
in (j) the relative viscosity may be expressed as 

T^r - V/Vs = 1 + 5 0h/2 + A{Pe, R/Rh) 0', (33) 

where 4>h = {Rh / FCf 4> is the volume fraction with 
respect to the hydrodynamic radius and the function 
A(Pe, R/Rh) contains the effects of microstructural dis- 
tortion (note that calculation of g(r) to 0(0) yields the 
viscosity to 0(0^)). The numerically determined func- 
tion A{Pe, R/ Rh) niay be further split into hydrody- 
namic and Brownian contributions, A = Ah + Ab, which 
are independently accessible from the numerical calcula- 



tions of [138|. 

The top panel of Fig[TT] shows A as a function of Pe. 
For shear rates up to Pe ~ 1 the qualitative variation 
of the viscosity is independent of the value of the size 
ratio R/Rh, displaying a low-shear Newtonian plateau 
followed by shear thinning. The R/Rh independence of 
the form of the curves for Pe < 1 reflects the fact that 
hydrodynamic interactions are not central to the mech- 
anisms underlying shear thinning and only influence the 
absolute value of the viscosity. For Pe>l the viscosity 
begins to increase as a function of Pe and the disper- 
sion shear thickens. In contrast to the shear thinning 
behaviour, the viscosity increase is strongly sensitive to 
the value of R/Rh- In the Limit R/Rh — >■ 1 the trajecto- 
ries sketched in FigHOlare recovered and Batchelor's pure 
hydrodynamic limit is realized with a viscosity indepen- 
dent of Pe. In |138| it is also shown that for R/Rh 1 
the normal stresses also vanish as Pe — > oo, indicating 
a Newtonian response. As R/Rh is increased the mag- 
nitude of the shear thickening reduces strongly and by 
R/ Rh = i-i is lost entirely. This trend strongly indicates 
the important influence of short range lubrication forces 
on shear thickening (see section [5?3| , which can effectively 
be turned-off by slightly reducing the hydrodynamic ra- 
dius below that of the repulsive potential interaction. 

The central and lower panels of Fig llll show the in- 
dividual hydrodynamic and Brownian pair contributions 
to the total stress as a function of Pe. For Pe < 1 it is 
apparent that the shear thinning is due to a reduction in 
the Brownian contribution with increasing Pe. For val- 
ues of R/Rh close to unity the reduction in As is more 
than compensated by an increase in the hydrodynamic 
stress afor Pe > 1, leading to shear thickening. However, 
increasing R/Rh above unity rapidly supresses the influ- 
ence of lubrication and the hydrodynamic contribution is 
overwhelmed by the strong drop in Ab- 

Comparing the flrst panel of FigHT] with FiglH it is 
remarkable the extent to which the qualitative rheologi- 
cal response observed in systems at finite volume fraction 
is reproduced by calculations based on the dilute limit. 
However, the fact that Wi = Pe at low volume fractions 
docs not permit investigation of potentially interesting 
interaction effects between shear thinning and thicken- 
ing. Despite the extensive understanding of the response 
of dilute dispersions to steady fiows, analogous solutions 
for time-dependent flows remain to be investigated. This 
leaves open many interesting questions regarding tran- 
sient response and non-steady states. 



6.3. Intermediate volume fraction 

The simplest way to extend the dilute-limit results to 
finite volume fraction is via the introduction of empiri- 
cal v olume fraction dependent scale factors [13, Il52l Il6ll 
Il67l |. Analysis of Eq.(|30l) has provided two important 
insights. Firstly, at finite volume fractions the value of 
g{r) outside the boundary layer close to the surface of 
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a reference particle should asymptote to the solution of 
V • [v(r) g{r)] = 0, for which the no-flux boundary con- 
dition is ignored. This amounts to assuming that, out- 
side the thin boundary layer, it is sufficient to solve a 
purely advective problem. Secondly, that the appropri- 
ate Peclet number in the presence of hydrodynamics is 
Pen = / 2D s{4>) , where Ds{(f>) is the short time dif- 
fusion coefficient (which differs from Dq due to hydro- 
dynamic interactions with neighbouring particles). For 
weak flows (Pe <C 1) these ideas are manifest in a modi- 
fied perturbation to the quiescent pair correlations 



.9(1-) = 9cq{r) 



1 



(34) 



The finite volume fraction equilibrium radial distribution 
function 5cq('') is an external input to the theory and 
can be calculated using either simulation or equilibrium 
integral equation theory [l^. The function f{r) is deter- 
mined by substitution of ()34p into the dilute limit equa- 
tion pop , or its generalization for non-hard-sphere poten- 
tials. Many-body effects are thus included via the equi- 
librium radial distribution function and the short-time 
diffusion coefficient entering Pen- It should be noted 
that Brady's approach assumes an input gcq(r) which di- 
verges at random close packing {(f) w 0.64). 

Despite considerable success, the scaling approach suf- 
fers from two significant drawbacks: (i) the mathematical 
structure of the nonequilibrium part of the theory is that 
of the dilute system. Thus, regardless of rescaling, this 
approach does not admit the occurrence of possible ad- 
ditional physical mechanisms which may only occur as 
a consequence of cooperative behaviour at finite volume 
fraction, (ii) the equilibrium microstructure 5cq('') is re- 
quired as an external input and does not emerge from an 
approximate treatment of many-body correlation effects 
within the theory. These issues may be addressed by 
approaches which aim to approximate the triplet distri- 
bution function entering the pair Smoluchowski equation 
via the effective force between a pair of particles 
(p5)) . In order to arrive at a closed theory it is neces- 
sary to relate g'^'^'' (ri, r2, ra), either explicitly or implic- 
itly (by approximating weighted integrals over the triplet 
distribution), to g(r) and the pair potential v{r) using an 
appropriate 'closure' hypothesis. 



6.3.1. Superposition approximation 

Guidance in developing an appropriate closure rela- 
tion to treat non-equilibrium states is provided by ex- 
perience from equilibrium liquid-state integral equation 
theory. One of the earliest approximation schemes aim- 
ing to arrive at a closed equation for the equilibrium 
pair correlations was developed by Born and Green [l68l | 
who employed the Kirkwood superposition approxima- 
tion g(3)(ri,r2,r3) sa g{ri2)g{r23)g{ri3) in combination 
with the second memb er of the exact Yvon-Born-Green 
heirarchy fEq [57|) p^ . Il69l | . Numerical solution of the 



resulting Born- Green equation yields acceptable results 
only for weak coupling. The superposition approxima- 
tion is asymptotically correct for large particle separa- 
tions but is poor when particles come close to contact, 
leading to a failure of the Born-Green equation at inter- 
mediate volume fractions (a breakdown which was erro- 
neously taken as an indicator for the first order freezing 
transition of hard-spheres). Subsequent attempts have 
aimed to systematically improve upon the superposition 
approximation b y inc luding additional Mayer cluster di- 
agrams (see e.g. [l70l |). 

One of the earliest attempts to close the non- 
equilibrium Ea.(|24p using the superposition approxima- 
tion, albeit i n th e absence of hydrodynamics, was made 
by Ohtsuki [l7l| (self diffu sion was addressed using an 
anologous approach in [l72j |). Numerical solution of the 
closed integro-differential equation resulting from this ap- 
proximation was performed for charged hard-spheres at 
intermediate volume fraction. Although the theoretical 
results for the zero-shear viscosity were found to be in 
reasonable agreement with those of experiment, the pair 
correlation g(r) was found to be in considerable error. 
These findings are supported by the work of Wagner and 
Russel who i nvestigated a similar superposition based ap- 
proach [M1\ . 

The observed discrepancies in g{r) arising from super- 
position are not surprising: In equilibrium, the pair flux 
in Eq. (j24[) may be set equal to zero, resulting in an exact 
equation for the pair correlation function 



fcsrV,lngeq(r) =F(r), 



(35) 



where F(r) is given by Eg. ([25]) and contains the unknown 
triplet distribution function. Use of superposition to ap- 
proximate ^('^^(ri, r2, ra) in Ea. (P5|) leads directly to the 
Born-Green equation for the equilibrium correlations, the 
shortcomings of which have been noted above. It is there- 
fore nontrivial that, despite a relatively poor description 
of the microstruct ure, t he results for the zero-shear vis- 
cosity presented in |l7l| turn out to be rather good agree- 
ment with experimental data. A similar situation is en- 
countered in a number of the theories to be outlined in 
this section and serves to highlight the fact that good 
values for integrated quantities (e.g. the viscosity) does 
not neccessarily imply that the underlying correlations 
are treated adequately. 

When applying superposition to tackle non- 
equilibrium problems it should be borne in mind 
that the approximation represents an uncontrolled 
ansatz and only possesses a firm statistical mechanical 
basis in the limit of vanishing flow rate. In equilibrium, 
the superposition approximation represents both the 
exact low volume fraction limit of the triplet correlation 
function and recovers correctly the long range asymp- 
totic behaviour. Analogous limiting results for the 
nonequilibrium triplet correlations which could motivate 
a more appropriate superposition-type approximation 
are currently lacking and would require a detailed anal- 
ysis of the triplet Smoluchowski equation in the dilute 
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limit. Despite these shortcomings, superposition may 
nevertheless turn out to be useful for some applications. 
As noted at the end of section [521 very recent simulation 
results investigating the microstructure of hard-sphere 
fluids under shear flow reve aled the importance of linear 
triplet configurations jl09l |. From equilibrium studies 
it is known that the error of the superposition approx- 
imation is reduced for such linear configurations |173| . 
thus raising the interesting possibility that superposition 
may be appropriate for treating certain non-equilibrium 
states. 



6.3.2. Potential of mean force 

Inspection of the exact expression for the force ((25|) 
shows that only a certain weighted integral of the triplet 
distribution is required to determine the pair correla- 
tions. It is thus not strictly neccessary to know the 
full details of the triplet distribution and schemes can 
be developed which aim to approximate directly the in- 
tegrated quantity. The simplest approximation possible 
is to neglect entirely the integral contibution to the ef- 
fective force (1251) and set 



F(r) = ~\Irv{r). 



(36) 



Combining this crude approximation with Eq. (j24p recov- 
ers the dilute limit equation of motion for g{r) for an 
arbitrary spherically symmetric pair potential v{r). For 
the special case of hard-spheres this leads to Eq. (pO]) , but 
with an additional delta function term on the right hand 
side. As demonstrated by Cichocki jl74l | the resulting 
equation is completely equivalent to (j30[) . with zero right 
hand side, supplemented by a no-fiux boundary condition 
at contact. 

As already noted, in equilibrium the pair Smolu- 
chowski equation ([24]) reduces to the Yvon-Born-Green 
equation for the pair correlations 



A:BTV^ln5eq(r) = -V,i;(r) 



x(V2w(|r2 
-Viw(r), 



2 7 goq(ri,r2) 
ral) - Viv(|ri -ral)) 
(37) 



where Wnif(r) is the equilibrium potential of mean force 
[iBl, defined by ^oq = exp(— /3wmf) in analogy with the 
low density limit of the pair correlations. A first step 
towards improving the zeroth order approximation (|36p 
is thus to approximate the nonequilibrium force psp by 
the equilibrium potential of mean force, leading to 



F(r) = -fcsrV,lngeq(r). 



(38) 



This approximation, developed by Russel and Gast 
[l48l |. incorporates equilibrium thermodynamic many- 
body couplings but, as is clear from Eq. ([57)) . neglects the 
influence of flow on the triplet correlations. It should be 



noted that the approximation (|38p does not provide any 
information regarding the nonequilibrium triplet correla- 
tion function. This is in contrast to superposition based 
approaches from which the triplet function can be recon- 
structed using a product of the self-consistently deter- 
mined pair correlation functions. It should also be noted 
that within the Russel-Gast approach the function gcqir) 
is an input, which can be calculated using either simula- 
tion or equilibrium statistical mechanical approximations 
(see subsection 16. 3. 3p . 

The linear equation resulting from combining equa- 
tions and ([55]) has been solved for hard-spheres 
in weak shear flow (l48j |. In these calculations, sim- 
ple approximations were employed to determine the hy- 
drodynamic functions B, G and H entering equations 
(|26| and (f28| defining the conditionally averaged hydro- 
dynamic tensors which represent the effective medium 
(solvent -I- {N — 2) colloids) in which the chosen pair of 
particles are immersed. Results were obtained for the 
zero-shear viscosity, linear response moduli G' (uj) , G" {uj) 
under small amplitude oscillatory shear and the leading 
order flow induced distortion of g{r) (via determination 
of the function /(r), see Eq. ipT]) ). For < 0.3 good 
agreement with experiment was obtained for the inte- 
grated quantities r]o,G' and G" , despite providing only 
a po or description of the microstructure (see Fig. 6 in 

An interesting feature of the Russel-Gast theory is that 
the predicted zero-shear viscosity becomes very large in 
the vicinity of random-close-packing {(j> « 0.64), although 
the precise nature of this rapid increase as a function of 
volume fraction remains to be studied in detail. The 
apparent divergence of 770 is a non-trivial output of the 
theory, given that the approximate Verlet- Weiss expres- 
sion for geq{r) used as input diverges only at (j) = 1.0 
[l76j |. When viewed within the context of the time- 
correlation/Green-Kubo formalism (see subsection 15. ip 
it is tempting to infer that the observed growth in 770 is 
related to the development of an underlying slow struc- 
tural relaxation time. The solution of Ea. (p8| for small 
amplitude oscillatory shear at finite frequencies would 
enable this issue to be addressed. 



6.3.3. Equilibrium integral equations 

The Russel-Gast theory outlined above neglects the in- 
fluence of external flow on the triplet correlation function, 
which leads to a force F(r) generated from the equilib- 
rium potential of mean force fEg. psp ). In order to go 
beyond this closc-to-equilibrium ansatz it is neccessary 
to express the force r(r) as a functional of g(r), such 
that both functions can be determined self-consistcntly. 
A promising approach in this direction is to generalize 
equilibrium liquid-state integral equation theory [l^ to 
treat the nonequilibrium problem. It will thus be use- 
ful to review briefly some concepts from the equilibrium 
theory before moving on to more unfamiliar territory in 
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the next subsection. 

Integral equation theories generally aim to calcu- 
late the equilibrium pair correlation function gcqir) 
from knowledge of the interaction potential in a non- 
perturbative fashion. Fundamental to the integral equa- 
tion approach is the Ornstein-Zernicke (OZ) equation 
which, for a translationally invariant system, is given by 
the convolution form Il6i 



(ri2) = Ccq(ri2) + n j drz Coq(ri3)/ioq(r32), (39) 

where hcq{r) ~ gcq{r)~l. The direct correlation function 
Ccq(r) defined by p9p is a function of simpler structure 
than heq{r) and is thus easier to approximate. When sup- 
plemented by an independent closure relation between 
Ceq(r) and hcq{r)^ containing details of the interaction po- 
tential under consideration, p9p provides a closed equa- 
tion for the pair correlations ( note that for pairwise addi- 
tive potentials the triplet distribution does not enter ex- 
plicitly) . The task of the integral equation practitioner is 
thus to find numerically tractable closures which capture 
the essential physics of the problem under consideration. 

Using diagrammatic techniques |l77| it can be shown 
that a formally exact closure relation is given by 

9cq{r) = exp[-/3w(r) + hcq{r) - Ccq(r) -f- bcq{r)], (40) 

where 6oq(r) is the unknown 'bridge function' containing 
the difficult to evaluate 'irreducible' Mayer cluster di- 
agrams [iBl- Two closures of particular merit are the 
Hyper- Netted-Chain (HNC) and Percus-Yevick (FY), 
given by 



bcq = 



bcq = -ln(l + /leq-Ccq) - (/l, 



cq '-•cq 



(HNC) (41) 
(FY) (42) 



respectively. An additional practical advantage of the- 
ories based on p9p over superposition-type approaches 
is that the convolution form of the integral term enables 
self consistent solutions to be obtained using efficient it- 
erative numerical algorithms jl78l | . 

Although the majority of integral equation theories fo- 
cus on the pair correlations, triplet corre latio ns can also 
be handled within the same framework |154| . In addi- 



tion to providing a higher level of resolution, the devel- 
opment of triplet-level integral equations is motivated by 
the desire for an improved description of the pair corre- 
lations. Exact relations, such as the YBG equation ([57|. 
connect the triplet to the pair correlations and the ex- 

(3) 

pectation is that errors in an approximate gcq may be 
averaged out by integration to the pair level. During 
the mid-1960s many of the leading liquid-state theorists 
proposed integral equations for the trip let correl ation s 
(e.g. Verlet UlSiM], Wertheim [Hi, Baxter [l8|, 
Stell |l86t ). all of which showed considerable promise. 
However, the complexity of solving the equations has 
hindered progress along this route and somewhat sim- 
pler, numerically tractable, theories now seem preferable 
I15l[l87l.[l88l. 



An integral equation which will be of particular rel- 
evance for t he following section was derived by Scher- 
winski [l89j |. Within the Scherwinski approximation 
the triplet correlation function is obtained from self- 
consistent solution of the following linear equation 

goq(ri,r2,r3) = gcq{r 12)9 cq{r 13)9 oq{r 23) +ngcq(ri2) X 
^ /^ .gcq (1-1,1-3, r4) , \\, / \l, f \ 

X / dr4 I ^ 5oq(^13) /lcq(?'14)/lcq(r24)- 



5cq(ri4) 



(43) 



Iteration of Eq. (|43| yields an infinite series expressing the 
triplet correlation function as a functional of the pair cor- 
relation function 5cq('')- Substitution of this series into 
the exact YBG equation ([57)1 yields a diagrammatic ex- 
pansion for 9cq{r) in perfect agreement with that arising 
from solution of Eqs.(13Hl), (001) and (gT]). In this sense, 
Ea. (|43p represents a triplet generalization of the more 
familiar pair-level HNC theory. The mo re powerf ul ap- 
proximations proposed by earlier workers |l79l - ll86 prob- 
ably provide a superior description of the triplet correla- 
tions than Eq. (|43p and would, upon substitution into the 
YBG equation, lead to improved (i.e. better than stan- 
dard HNC) estimate of the pair correlations. However, 
the Scherwinski approximation has a number of purely 
technical advantages which make it particularly suitable 
for application to non-equilibrium situations and which 
are convenient for nu meri cal implementation (see also 
jl47l |. which predates [l89j |. but contains several of the 
key ideas). 



6.3.4- Noneqmhbrium integral equations 

In order to go beyond the Russel-Gast approximation 
jl48j outlined in subsection I6.3.2[ Lionberger and Russel 
employed the Scherwinski equation for the triplet corre- 
lation function (|43| in order to estimate the force (|25p en- 
tering the pair Smoluchowski equation [l90j j . Nonequilib- 
rium pair and triplet correlations are thus determined self 
consistently and are both influenced by the externally im- 
posed flow. It should be noted that the direct application 
of an equilibrium relation, such as Eq. pS)) . to nonequi- 
librium ignores the fact that nonequilibrium states are 
intrinsically different from equilibrium and thus repre- 
sents a major approximation. 

The origin al ve rsion of the Lionberger-Russel theory 
presented in [l90j | neglected hydrodynamic interactions 
and has been implemented numerically for weak flows 
only. No results beyond leading order in Pe have been 
presented, although in principle the theory remains valid 
also in the nonlinear regime. For < 0.4 the LR the- 
ory makes predictions for the zero-shear viscosity, self- 
diffusion cooefhcient and distorted microstructure in rea- 
sonable agreement with available computer simulation 
results. For (f> > 0.45 significant quantitative deviations 
appear and the theory becomes unreliable. We note that 
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the study [190l | considered a suspension interacting via 
a continuous repulsive potential which was then mapped 
onto a hard-sphere system. The input equilibrium mi- 
crostructure wa s gen erated using the Rogers- Young in- 
tegral equation [191| (an interpolation between PY and 
HNC), despite the fact that the Pe — > limit of the 
theory reduces to the HNC approximation for goq{f)- 

Although the Lionberger-Russel theory |l9C| provides 
a sophisticated treatment of the microstructural distor- 
tion by incorporating the triplet correlations into the 
self-consistency loop, it is interesting that the results for 
rjo are inferio r to those from the much simpler Russel- 
Gast theory |l48| at high volume fractions. In partic- 
ular, the former predicts only a relatively weak growth 
of rjQ with volume fraction, whereas the latter suggests a 
divergence. Calculations of the linear moduli G',G" as 
a function of frequency have been performed using the 
LR theory and, perplexingly, reveal a structural relax- 
ation time which decreases with increasing volume frac- 
tion. This unphysical prediction would appear to be at 
odds with the growth of 770 output from the theory and 
represent s a w eak point of the approach, in need of clar- 
ification |l92| . We note that the only approximation in- 
voked by the Lionberger-Russel theory is the Scherwinski 
closure for the triplet correlation function. It would 
therefore be interesting, albeit numerically demanding, 
to see whether any of th e more sophisticated triplet clo- 
sures available [T79l - [l88t can improve the performance at 
higher volume fractions. 

The theory developed in |l90l | omitted hydrodynamic 
interactions in order to make clearer the approximations 
to the many-body thermodynamic couplings and to fa- 
cilitate cor npar ision with Brownian dynamics simulation 
results. In [192l | hy drody namic interactions were included 
into the theory of jl9d | in order to make closer contact 
with experiment. The hydrodynamic approximations de- 
veloped were found to be of the correct magnitude but 
errors in the underlying thermodynamic approximation, 
namely the Scherwinski equation (|43p . led to an under- 
estimation of the magnitude of the nonequilibrium struc- 
ture. It was identified that the magnitude of the flow 
induced structural distortion is determined by the slow 
structural relaxation in the syste m (a finding supported 
by subsequent theoretical studies [l93j), which, taken to- 
gether with the results for the volume fraction depen- 
dence of the linear response moduli, suggests that the 
integral equation approach docs not capture the slow dy- 
namics characteristic of dense systems. 



Wagner and Russel |147l | developed an integral equa- 
tion approach closely related to that of Lionberger and 
Russel in which the nonequilibrium triplet correlations 
are approximated using a closure motivated by the PY 
equilibrium theory P^ . Although use of a PY-type clo- 
sure ensures that excluded-volume packing constraints 
are treate d rea listically, at least close to equilibrium, the 
theory of [147{ included hydrodynamic interactions by 
employing only the low density limit of the pair Hydro- 
dynamic functions A, B,G and H entering Eqs. (|26p and 



(p8|) . The simultaneous introduction of hyd rodynamic 
and thermodynamic approximations in jl47| served to 
obscure the validity of the proposed nonequilibrium PY 
approximation. 

When using integral equation methods to tackle the 
triplet corclations in nonequilibrium it is important to 
bear in mind that the physical situation is intrinsically 
different from that in equilibrium. Consequently, cau- 
tion must be exercised when attempting to apply trusted 
and familiar results from equilibrium statistical mechan- 
ics to system under flow. In order to appreciate more 
clearly the approximations involved in applying equilib- 
rium triplet closures to the pair Smoluchowsk i equ ation 
both Lionberger and Russel |l90l | and Szamel |175| have 
invoked the concept of a 'fictitious' flow dependent two- 
body potential w(r,7). In a study of the k inetic theory 
of hard-spheres Rcsibois and Lcbowitz [l94l | assumed the 
existence of a two-body potential u(r,7) which, if em- 
ployed in an equilibrium calculation, reproduces exactly 
the nonequilibrium pair correlation function 



= ,9cq(r; [u]), 



(44) 



where the square brackets indicate a functional depen- 
dence and where geq is anisotropic as a result of the aniso- 
topy of u. The fictitious potential thus serves as proxy 
for the flow field acting on the real system. It should 
be made clear that the assumption that an effective two- 
body potential can yield the correct g{r) is quite distinct 
from the (erroneous) assumption that an effective one- 
body external potential field can represent the one-body 
flow induced force acting on the particles (see the end of 
section [5?^ . Eq.(|33]) implicitly assumes that a homoge- 
neous one-body density distribution n is not a function 
of the flow rate, thus neglecting possible dilation effects. 
Given equation it is natural to go one step further 
and assume that the relation can be uniquely inverted, 
such that 



u(r,7) = u(r; [g]). 



(45) 



By construction, u(r, 7) reproduces the nonequilibrium 
pair correlations. However, if the same fictitious two- 
body potential is used in an equilibrium statistical me- 
chanical calculation of the triplet correlation function, 
the exact nonequilibrium g^"^^ will not be reproduced. 
The missing part of the triplet correlation is referred to 
as the 'irreducible' term 

(r 1 , ra , rs , 7) = 5^^) (^ , rs , rg ; [ 5 ] ) + g^^^^ (r 1 , ra , rg ; 7) , 

(46) 

where we have assumed that l|45p is valid. A key approx- 



,(3) 



0, 



imation in the work of |19G | and jl7a | is to set gl 
which essentially amounts to assuming that equilibrium 
relations such as (|43)) may be used to connect triplet and 
pair functions in nonequilibrium. 

The nonequilibrium integral equation method consid- 
represents a synthesis of the exact 
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dilute-limit results with statistical mechanical descrip- 
tions of dense systems at equilibrium. While this ap- 
proach is promising for weak flows {Pe <C 1), data is 
lacking for stronger flows which would enable the close- 
to-equilibrium character of the closure approximations to 
be better tested. In particular, it has been emphasized 
by Szamel that a major challenge for theories of 

the distorted microstructure is to account for the nonzero 
distortion of ^(r) in the vorticity-gradient plane perpen- 
dicular to an applied shear flow. This distortion cannot 
be detected to leading order in Pe and makes desireable 
numerical studies of the existing closure approximations 
under strong flow. Applications of the nonequilibrium 
integral equation method have been restricted to shear 
flow and, wi th the exc eption of small amplitude oscil- 
latory shear |l90l . Il92j . have neglected time-dependent 
flows entirely. 

6.3.5. Alternative approaches 



1 
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FIG. 12: The reduced zero-shear viscosity as a function of 
volume fraction as predicted by various th eories based on the 
pair Smoluchoski equation : Sz amel (black) [175||, Br ady (red), 
Lionberger-Russel (blue) [190| |. Ronis (green) [l95j . The in- 
set shows the results of the Szamel theory compared with 
Brownian dynamics simulation data (circles). For (j) > 0.43 
the theory strongly und erestimates the zero-shear viscosity. 
(Adapted from Ref.[l7i|) 



The nonequilibrium pair correlations and rheology of 
colloidal dispersi ons under weak shear flow was investi- 
gated by Szamel |175j , who employed functional methods 
to approximate the unknown integral term in Ea. (|25p . 
Hydrodynamic interactions were neglected entirely. As- 
suming that the irreducible term in (pS)) can be neglected, 
the triplet correlation function can be developed in a 
functional Taylor expansion about the equilibrium state. 
To flrst order this is given by 



g(3)(r^^r2,r3;[5]) =5^q^(ri,i-2,r3) 



(3)/ 



(47) 



dYi dr5 — — ^ — g i-4,r5 - ffoq , 



where the functional derivative is evaluated for a con- 
stant one-body density. Equation ()47|) is clearly a close- 
to-equilibrium approximation. Insertion of (|47p into the 
expression for the force (j25p and rearrangement of terms 
leads to a closed equation for ^(r) which requires both 

ffoq (ri2,^23,?'i3) and c^qVi, r2, i"3, r4), a higher order 
equilibrium direct correlation function, as input. Fol- 
lowing appropriate decoupling approximations for the 
unknown higher order equilibrium correlations and lin- 
earizing with respect to the shear flow, Szamel obtained 
a closed equation for ^(r) which only requires gcqir) as 
input ( for w hich the Verlet- Weiss approximation was em- 
ployed [nl). 

The Szamel theory [l75j is considerably simpler to im- 
plement than the integ ral equation approach of Lion- 
berger and Russel |l90l | and provides comparable, per- 
haps even slightly better, results for the zero-shear vis- 
cosity as a function of volume fraction, as is evident 
from FiglT^l Both g{r) and its Fourier transform ^(k) 
were found to be qualitatively similar to those from the 
Lionberger-Russel theory, substantially underestimating 
the magnitude of the distortion from equilibrium when 
compared to Brownian dynamics simulation. 



In an early study, Ronis took an alternative ap- 
proach based on fluctuating hydrodynamics, in which 
phcnomenological fluctuating terms are add ed to the 
macroscopic equations of hydrodynamics |195j . The the- 



ory of stochastic processes may then be employed to cal- 
culate nonequilibrium time-correlation functions and the 
distorted microstucture. Of all the theories described in 
this and the previous section, Ronis theory provides the 
most accurate results for the distortion of the structure 
factor at low Pe values. Nevertheless, it can be shown 
that the Ronis approximation leads to a vanishing of the 
microstructural distortion in the vorticity-gradient plane 
at all Pe, i n co ntradiction to experiment and simula- 
tion results jl66j . This deflcit is related to the fact that 
the hard-sphere 'core condition' g{\r\ < 1) = is vio- 
lated within this approach. The Ro nis th eory reduces to 
the closely related theory of Dhont jl63l | in the low Pe, 
linear ized limit. In another early work, Schwarzl and 
Hess |l96| postulated a phenomenological equation for 
g{r) involving a number of empirical parameters repre- 
senting the relaxation times in the system. However, due 
to the phenomenological nature of both the fluctuating 
hydrodynamics approach and the equation-of-motion for 
g(r) proposed by Schwarzl and Hess, t he f ound ation of 
the theoretical treatments presented in |195| and [l96j in 
nonequilibrium statistical mechanics is unclear. 

Finally, we note that Wagner has assessed approaches 
based on the pair Smoluchowski equation using the 
GENE RIC f ramework of beyond equilibrium thermody- 
namics jl97l |. This formalism enables any proposed clo- 
sure of the pair Smoluchowski equation to be checked 
for the rmodynamic consistency. The study presented 
in jl97| identified the thermodynamically admissable ex- 
pression for the stress tensor and clarified the nature of 
the inconsistencies which can occur when separate deriva- 
tions arc performed for the equation-of-motion for g{r) 
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and the stress tensor. 



6.4. Temporal locality vs. memory functions 

In sections f5.1l and [5T2] wc introduced briefly the Green- 
Kubo expressions for calculating the shear stress in both 
the hnear ([T7| and nonlinear regimes ([22]) . Within the 
Green-Kubo framework, transport coefficients are related 
to integrals over time correlation functions jl99l |. The 
growth of the shear viscosity as a function of volume 
fraction is thus related to an increasingly slow decay 
of the stress autocorrelation function. What is lack- 
ing thus far in the discussion of the pair Smoluchowski 
equation and its approximate solutions, is the connec- 
tion between the temporally local Eas. ([M|) and ([25)1 and 
the nonlocal expressions ([TT]) and (|22p fundamental to 
the time-correlation function formalism. Moreover, the 
memory kernels characteristic of the time-correlation ap- 
proach (and widely employed in continuum mechanics 
approaches, see section [2]) make no explicit appearance 
within the pair Smoluchowski framework. 

The key to understanding the connection between 
the nonlocal time-correlation formalism and approaches 
based on the local pair Smoluchowski equation lies in 
the study of time-dependent external flow fields. In- 
deed, it is the absence of time-dependent data from the 
pair Smoluchowski approach which has served to ob- 
scure the relationship between these two methods, de- 
spite the fact that they are formally equivalent. The 
only available Smoluchowski-based time-dependent cal- 
culations were perf ormed using the Lionberger-Russel 
theory [281 [iM Il92| for small amplitude oscillatory shear 
flow. As pointed out in subsection 16. 3. 41 the available re- 
sults for the volume fraction dependence of G'{uj) and 
G"(u)) within the LR theory reveal underlying problems 
resulting from the approximations employed (i.e. a re- 
duction of Tq, with increasing 0) which would otherwise 
have gone unnoticed. This serves to emphasize the im- 
portance of going beyond steady flow calculations in ap- 
plications of pair Smoluchowski theories. 

The generalization of Eq. ip^ to general time- 
dependent shear is given by integration over the entire 
flow history [l^ 

axy(t) = / dt'jit')G{t,t'), (48) 

^ — oo 

where G{t,t') = G{t,t'; [7]) is the nonlinear shear mod- 
ulus. The lack of time-translational invariance in the 
modulus arises from a functional dependence on the shear 
rate. The microscopically derived Eq. p8|) should be con- 
trasted with the more familiar phenomenological result 

Using the generalized Green-Kubo result (|48| it is in- 
structive to consider a simple special case: The stress re- 
sponse of hard-spheres to the onset of steady shear flow 
in the absence of hydrodynamic interactions. Specifically, 



we consider a shear field which is switched from zero to 
a constant value, j{t) = 70(— <). For this choice of shear 
field Eq.dH]) reduces to 

<J^y{t)=i f dt'Gssit'), (49) 
Jo 

where Gss{t) is the time translationally invariant shear 
modulus under steady shear flow. On the other hand, 
within approaches based on the distorted microstructure 
the interaction contribution to the stress is given by (see 
Eq.(l29l)) 

Ti^ f TV 

^^yit) = - — jdv-v\r)g{v,t), (50) 

where g{v,t) is calculated from the pair Smoluchowski 
equation ([24| subject to the switch-on shear flow under 
consideration. The expressions (|49l) and (|50l) are formally 
equivalent. Equating the time derivatives thus leads to 
the exact relation 

Gssit)=-—Jdr-v{r)^^. (51) 

The quiescent shear modulus is thus recovered in the slow 
flow limit 

The right hand side of Eq. ([5^ can be further reduced to 
the determination of the function f(r,t) by substitution 
of the first order expansion (|3ip for g{r,t). 

As the volume fraction is increased, the time deriva- 
tive on the r.h.s. of Eg. ((5^ must give rise to a growth in 
the timescale determining the relaxation of stress fluctu- 
ations, as described by Geq(t). The fact that this is not 
captured by the Lionberger-Russel theory [H, Il90l . Il92j 
simply reflects the failings of the approximate closure 
(|43p relating the nonequilibrium triplet and pair corre- 
lation functions. However, it is not clear that even an 
exact equilibrium expression for the triplet correlations 
would be sufficient to resolve these difficulties. If the key 
source of error lies in the neglect of the irreducible term 
apprearing in (|46|) then considerable new insight into the 
nature of nonequilibrium states will be required in order 
to make further progress using integral equation meth- 
ods. Nevertheless, these considerations should ser ve to 
motivate ti me-de pendent studies of the Russel-Gast |148| 
and Szamel |l75l | theories. 

While the particular example chosen (switch-on shear 
fiow) provides access to the steady shear rate dependent 
modulus, other choices of time-dependence may enable 
connection to be made between non-time translationally 
invariant correlation functions (e.g. G{t,t')) and time- 
dependent solutions of the pair Smoluchowski equation. 
The issue of whether temporally local constitutive equa- 
tions are preferable to nonlocal functionals for describing 
complex fluids has been addressed within the framework 
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of nonequilibrium thermodynamics |198j . The nonequi- 
Hbrium thermodynamics approach requires identification 
of an appropriate set of structural variables which, in ad- 
dition to the standard hydrodynamic variables of mass, 
momentum and internal energy density, contain all in- 
formation about the state of the system at a given time 
that is necessary to determine the macroscopic quantities 
of interest. The pair Smoluchowski approaches discussed 
in this section essentially introduce gir) as an additional 
structural variable \l9l\ . While this appears to be a valid 
approach for low and intermediate volume fractions, a 
correct identification of the structural variables appropri- 
ate for describing glass formation and dynamical arrest 
remain to be found. For this reason the most promis- 
ing approaches to treating high volume fraction states 
arc based on the generalized Green-Kubo relations and 
mode-coupling theory. 



7. GLASS RHEOLOGY 

Assuming that crystallization effects can be suppressed 
(see section U) the volume fraction can be increased to 
the point at which the individual particles are unable 
to diffuse beyond the cage of nearest neighbours and 
a dynamically arrested glassy state is formed. In or- 
der to visualize the amorphous cage structure in such 
a glassy state Figure [13] shows a configuration snapshot 
taken from a Brownian dynamics simulation of a bi nary 
hard-sphere mixture in two dimensions (hard-discs) [T^l ■ 
The two-dimensional volume fraction (/)2d = 0.81 of the 
simulation is above the estimated glass transition point 
of 02D ~ 0.79 and the size ratio of large to small disc 
radii is 1.4, a value empirically found to frustrate crys- 
tallization in two dimensional systems (which occurs at 
'?^2D = 0.69 for monodispcrse discs). In both two- and 
three-dimensional systems the physics of the glass tran- 
sition becomes important for determination of both the 
rheology and flow distorted microstructure of high vol- 
ume fraction systems. 

The response to externally applied flow of states close 
to, or beyond, the glass transition is only beginning to be 
understood and establishing the basic principles of glass 
rheology remains a challenging task. At present, the only 
truly microscopic theories available are provided by re- 
cent extensions of the quiescent MCT which enable the 
effects of external flow to be incorporated into the formal- 
ism and thus make possible a theoretical investigation of 
the complex interaction between arrest and flow. 



7.1. MCT inspired approaches 

Extending earlier work on the low volume fraction self 
diffusion of colloidal dispersions |200l | , Miyazaki and Re- 
ichman constructed a self-consistent mode-coupling-type 
approach to describe collective density fluctuations for 
dense colloidal fluids under shear below the glass transi- 




FIG. 13: A snapshot from a Brownian dynamics simulation of 
a quiescent binary hard-disc mixture (using a size ratio of 1 : 
1.4 to supress crystallization). The simulation was performed 
at a two-dimensional volume fraction of (j>2D ~ 0.81, which 



lies above the estimated glass transition packing 



a(9) 



' 0.79, 



with 50% large discs and 50% small discs. (Figure courtesy 
of F. Weyfier.) 



tion |20ll4203| . The Miyazaki- Reichman theory consid- 
ers time-dependent fluctuations about the steady state 
and thus requires the (unknown) flow distorted structure 
factor S'(k) as an input quantity. Approximating S'(k) 
by the quiescent correlator, results have been presented 
for co Uoidal dis persions in two-dimensions under steady 
shear [20lLl202l | and in three dimensions (subject to addi- 
tional isotropic approximations) under oscillatory shear 
[203j | . Applications to glassy states have been avoided as 
the theory relies upon an ergodic fluctuation-dissipation 
theorem. An alternative extended-MC T ap proach has 
been prop osed by Koblev and Schweizer 204l | and Saltz- 
man et al [205j | which is built on the idea that cntropic 
barrier hopping is the key physical process driving the 
microscopic dynamics and rheology of glassy colloidal 
suspensions. Due to the activated nature of the barrier 
hopping process the ideal glass transition described by 
quiescent MCT (see section|3|) plays no role. A nonlinear 
rheological response results from a stress induced modi- 
fication of the barrier heights. 

A currently promising method of extending quiescent 
MCT to treat dense systems under flow involves integra- 
tion through the transient dynamics, starting from an 
equilibrium Boltzmann distribution in the infinite past. 
In contrast to ^201, 203J the distorted microstructure is 
an output of this approach. The initial form of the thc- 
ory w as outlined by Fuchs and Cates for steady shear flow 
[l40l | and presented two essential developments: Firstly, 
that integration through the transient dynamics leads di- 
rectly to exact generalized Green-Kubo formulas, relat- 
ing average quantities to integrals over microscopic time- 
correlation functions. Secondly, that MCT-type projec- 
tion operator approximations reduce the formal Green- 
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Kubo expressions to closed equations involving transient 
correlators (which can be calculated self-consistently). A 
strong prediction of the ITT-MCT theory resulting from 
combining these two steps is that the macroscopic flow 
curves exhibit a dynamic yield stress (see section [5^ in 
the limit Pe— >0, for states which would be glasses or gels 
in the absence of flow. Moreover, the yield stress appears 
discontinuously as a function of volume fraction, in con- 
trast to mesoscopic approaches fsi-fSTt- The ITT-MCT 
thus provides a scenario for a nonequilibrium transition 
between a shear-thinning fluid and a yielding amorphous 
solid which is supported b y c onsi derable ev idence from 
both colloidal experiment s Ipp. \q^ . Il39l |206| and Brown- 
ian dynamics simulation j4lL l79l|. 

Due to the numerical intractability of the microscopic 
theory of |l4d | , subsequent work focused on the construc- 
tion of both isotropically averaged approximations to the 
full anisotropic equations and simplified schematic mod- 
els inspired by these [s^. Comparison of the theoret- 
ical predictions with experimental data for thermosen- 
sitive core-she ll particles ( see Figl?]) has proved highly 
successful [63l.l64. [207l - [209j . The original formulation of 
the ITT-MCT (more details of which can be found in 
[210l |') has subsequently been superseded by a more ele- 
gant version [2l| . It is interesting to note that the signif- 
icant technical changes to the ITT-MCT formalism in- 
troduced in pTj lead to expressions whi ch resemb le more 
closely those of Miyazaki and Rcichman |20l| - |203| . Given 
the very different nature of the formal derivation (fluctu- 
ating hydrodynamics vs. projection operator methods) 
and approximations employed, the similarity of the fi- 
nal expressions is reassuring and serves to highlight the 
robustness of MCT-based approaches. For a comprehen- 
sive overview of the status of the steady shear theory we 
refer the reader to the recent review |20j . 

Going be yond steady shear, the original formulation of 
ITT-MCT [HiliiSIilfll has been generalized to treat ar- 
bitrary time-dependent shear \vi\ . These developments 
not only enable shear fields of particular experimental 
relevance to be investigated (e.g. large amplitude oscil- 
latory shear flow), but have also revealed an underlying 
mathematical structure which is not apparent from con- 
sideration of steady flows alone. The theory has been 
applied (albeit subject to various simplifying approxi- 
mations) to investigate the build-up of stress and corre- 
sponding microscopic particle motion, as encoded in the 
mean-squared-displacement, following the onset of shear 
[esf . More recently, the modern version of ITT-MCT [21| 
has been extended to describe time-dependent flow of ar- 
bitrary geometry [isj , thus making possible the study of 
non-shear flow and enabling the full tensorial structure of 
the theory to be identified. The developments presented 
in [ll] elevate the ITT-MCT approach to the level of a 
full constitutive theory for dispersion rheology and may 
be regarded as the most up-to-date formulation of the 
theory. While the development of numerical algorithms 
to efficiently solve the fully microscopic theory [l^ is 
currently in progress, this task is made computationally 



demanding by the combination of spatial anisotropy and 
lack of time-translational invariance presented by many 
flows of interest. In [l^ a simplified theory was pre- 
sented which contains the essential physics of the full 
microscopic equations, including the tensorial structure, 
but which is much more convenient for numerical solution 
(see subsection 17. 7p . 



7.2. Integration through transients 

The integration through transients ( ITT ) approach 
originally developed by Fuchs and Gates jl4d j and subse- 
quently developed in [ill provides a formal expres- 
sion for the nonequilibrium distribution function ^'(t) 
required to calculate average quantities under flow. In 
essence, ITT provides a very elegant method of deriving 
generalized (i.e. nonlinear in K(t)) Green-Kubo relations 
which invite mode-coupling-type closure approximations. 
The current formulation of the theory neglects hydrody- 
namic interactions (HI) entirely. On one hand this omis- 
sion is made for purely technical reasons, but it is also 
hoped that HI will prove unimportant for the microscopic 
dynamics of the dense states to which the final theory 
will be applied. In the following we will briefly outline 
the key points of ITT, employing throughout the modern 
formulation of [ill ■ 

The starting point for ITT is to re-express the Smolu- 
chowski equation (|14p in the form 



dt 



(53) 



where, in the absence of HI, the Smoluchowski operator 
controlling the dynamical evolution of the system is given 

by [a 



n{t) = ^ a, • [I?o(9, - /3F0 - K{t) • r,]. 



(54) 



Equation (j53p may be formally solved using a time- 
ordered exponential function (which arises beca use f l(t) 
does not commute with itself for different times |21l| ) 



cxpi 



ds il(s) 



(55) 



where 'I'oq is the initial distribution function in the 
inflnite past, which is taken to be the equilibrium 
Boltzmann-Gibbs distribution corresponding to the ther- 
modynamic state-point under consideration. The as- 
sumption of an equilibrium distribution is clearly accept- 
able in situations for which the quiescent state is one of 
ergodic equilibrium. However, the role of the initial state 
is less clear for statepoints in the glass and the depen- 
dence, if any, on the initial condition may depend upon 
the details of the flow between t = —oo and the present 
time t. The absence of a general proof that ^'(i) is inde- 
pendent of 4'(— cxd) leaves open the possibility that cer- 
tain flow histories do not restore ergodicity and that the 
system thus retains a dependence on the initial state. 
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The solution (|55|) is formally correct, but not partic- 
ularly useful in its present form. A partial integration 
yields an alternative solution of (|53|) which is exactly 
equivalent to (|55p. but more suited to approximation 

xSj{t) = / dti *eK;(ii):6-e_i , (56) 

^ — oo 

where = — J2i ^"''f ^'^'^ the 'double dot' notation 
fami liar from continuum mechanics, A : B = Tr{A • B} 
j212l |. has been employed. As a result of the partial in- 
tegration the dynamical evolution in Ea. (j56p is dictated 
by the adjoint Smoluchowski operator 

n\t) = J2[Do{d, +f3F,) + r, ■ K^{t)]- 8,. (57) 

i 

Equation ([5S)) is the fundamental formula of the ITT 
approach and expresses the nonequilibrium distribution 
function as an integral over the entire transient flow his- 
tory. Both solutions (|55p and ([55]) are valid for arbitrary 
flow geometries and time-dependence. The relation be- 
tween the two formal solutions is analogous to the Heisen- 
berg and Schrodinger pictures of quantum mechanics in 
which the time evolution of the system is attributed to 
either the wavefunction (equation (|55p ) or the opera- 
tors representing physical observables (equation ((56| V It 
should be understood that the ITT form ([55]) is an oper- 
ator expression to be used with the understanding that 
a function to be averaged is placed on the right of the 
operators and integrated over the particle coordinates. 
A technical point to note is that in cases for which phase 
space decomposes into disoint pockets ('nonmixing' dy- 
namics) the distribution (|56p averages over all compart- 
ments. A general function of the phase-space coordinates 
/(t, {r^}) thus has the average 

{fr = {f)+ / dh{K{U):ae^l^ V), (58) 

where (/)"'^ denotes an average over the non-equilibrium 
distribution ([56]). Equat ion ([ 56 ]) generalizes the orig- 
inal formulation of ITT pot to treat arbitrary time- 
dependence [isj . 

7.3. Translational invariance 

Before applying mode-coupling-type approximations 
to the exact result (|58p we first address an important con- 
sequence of assuming homogeneous flow (refiected in the 
spatial constancy nit) appearing in Eg. ilSSp ). On purely 
physical grounds, it seems reasonable that for an infinite 
system the assumed translational invarance of the equi- 
librium state (crystallization is neglected) will be pre- 
served by the Smoluchowski dynamics. However, prov- 
ing this for a general time-dependent flow is mathemat- 
ically nontrivial, due to the fact that the Smoluchowski 
operator ((M)) is itself not translationally invariant. By 



considering a constant vectorial shift of all particle co- 
ordinates, y[ = Vi + a, Brader et al have shown that 
the nonequilibrium distribution function Vl'(t) is transla- 
tionally invariant (but anisotropic) for any homogeneous 
velocity gradient Hi{t) (isj . 

Given the translational invariance of ^'(t) it becomes 
possible to investigate the more interesting invariance 
properties of the two-time correlation functions. The 
correlation between two arbitrary wavevector-dependent 
fiuctuations 5fc^ = /q - (/q)"'' and 5g\, = .gk ~ (ffk)"'' 
occuring at times t and t' is given by 

Cf^sStX)^{5r^{t)5g^{t')r. (59) 

It is clear that in a homogeneous system the correlation 
function (|59p must be translationally invariant. How- 
ever, in this case, shifting the particle coordinates by a 
constant vector a yields 

Cf^aSt^t') = e-'('»(*'*')-'')-C/^,Jt,i'), (60) 

where 

q(f,t') = q-e:^*''''"^'^ (61) 

The only way in which the required translational invari- 
ance of the correlation function can be preserved is if the 
exponential factor in (j60p is equal to unity. This require- 
ment has the consequence that a fluctuation at wavevec- 
tor k = q(t, t') at time t' is correlated with a fluctuation 
with wavevector q at time t as a result of the affine sol- 
vent flow. Ea. (|6ip thus deflnes the advected wavevector 
which is central to the ITT-MCT approach and which 
captures the afRne evolution of the system in approaches 
focused on Fourier components of fluctuating quantities 
(e.g. the density p^J rather than particle coordinates di- 
rectly. The wavevector q(t,<') at time t' evolves due to 
flow-induced advection to become q at later time t. It 
should be noted that various definitions and notations for 
the advected wavevector have been employed in the lit- 
erature documenting the development of ITT-MCT and 
which could provide a source of confusion. In the present 
work we exclusively employ the modern definition used 

in HIiiiil. 

Although Eq. (pT|) arises from microscopic considera- 
tions it is nevertheless fully consistent with the contin- 
uum mechanics approaches outlined in section [51 de- 
spite the very different mindset underlying the two meth- 
ods. Ea. (|5T|) simply describes the affine deformation of 
material lines in Fourier space and can thus be used 
to define the inverse deformation gradient tensor via 
q(t,t') = q • 'E~'^{t,t') in complete accord with contin- 
uum approaches. Doing so leads to the identification 

E-i(i,t')-e:^*"''"^'^ (62) 

As the deformation gradient tensor simply describes the 
affine distortion of a material line under flow, it is nat- 
ural to deflne also a reverse- advected wavevector result- 
ing from the inverse transformation q(t, t') = q • E(t, t'), 
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where 



E(t,t') 



J*,dsn{s) 



(63) 



The choice of using either advected or reverse-advected 
wavevectors in treating the effects of affine motion within 
a microscopic theory has parallels with the choice be- 
tween Lagrangian and Eulcrian specifications of the 
flow field in continuum fluid dynamics approaches [3l| . 
Within a continuum mechanics framework the deforma- 
tion gradient would simply be defined as the solution of 
the equation 



d_ 
dt 



E{t,t') = K{t)E{t,t'), 



(64) 



for a given flow K{t). Acco rding to the rules of time- 
ordered exponential algebra |21l| . Ea. (p5|) is the formal 
solution of (|64)). thus demonstrating the consistency be- 
tween the Fourier-space microscopic approach of |18| and 
traditional real-space continuum mechanics. 

The advected wavevector introduced above provides a 
convenient way to keep track of the affine deformation 
in a particulate system. Mode-coupling-type approxima- 
tions (to be discussed below) seek to factorize the av- 
erage entering Eq.(|58|) by projecting t he dy namics onto 
the subspace of density fluctuations Pq jl03| . For a flow- 
ing system a fluctuation at wavevector q{t,t') at time t' 
evolves (in the absence of interactions and Brownian mo- 
tion) to one at q at time t. It thus becomes essential to 
project onto density fluctuations at the correct advected 
wavevectors in order to avoid spurious decorrelation ef- 
fects in the resulting approximations. 



7.4. Microscopic constitutive equation 

In order to address dispersion rheology the special 
choice / = &/V is made in ([55]) . leading to an exact 
generalized Green-Kubo relation for the time-dependent 
shear stress tensor [la 



(T(t) = — / dti{K{ti):(T ej cr), 

*^ J — oo 



(65) 



noting that there are no 'frozen in' stresses in the equilib- 
rium state, (cr) = 0. The adjoint Smoluchowski operator 
^'^{t) has a linear dependence on K{t) and equation ((55)) 
is thus nonlinear in the velocity gradient tensor. Equa- 
tion (|65p is a formal constitutive equation expressing the 
stress tensor at the present time as a nonlinear func- 
tional of the flow history. Although the result (|65|) docs 
not provide an exact description of the physical system 
under consideration (particle momenta are assumed to 
have relaxed and hydrodynamic interactions are absent), 
it has a formal status equivalent to that of Eg . ([55)) . For 
the special case of steady shear flow is consistent 
with (22) with a shear modulus given by (P5)) . 



Application of MCT-type projection operator factor- 
izations [3 to the average in (p5)) leads to a compli- 
cated, but closed, constitutive equation expressing the 
deviatoric stress in terms of the strain history [l^ 



a{t) 



dt' 



dk 
32^ 



-(k-B(i,i')-k)kk 



kk{tX}Sl 



k(t,t') 



{t,t') 



(66) 



where Sk and 5^ are the equilibrium static structure 
factor and its derivative, respectively. The influence of 
external flow enters both explicitly, via the Finger ten- 
sor B(t,i') (see subsection 12. ip . and implicitly through 
the reverse-advected wavevector. As noted above, the 
reverse-advected wavevector, which provides an impor- 
tant source of nonlinearity in (|66p . enters as a result of 
judicious projection of the dynamics onto appropriately 
advected density fluctuations /Ok(t,t')- The normalized 
transient density correlator describes the decay under 
flow of thermal density fluctuations and is defined by 



*k(i,i') 



(67) 



The occurance of the advected wavevector in (|67|) en- 
sures that trivial decorrelation effects are removed (i.e. 
that in the absence of Brownian motion and potential 
interactions $k = 1 for all times). 

In order to close the constitutive equation (|66| we re- 
quire an explicit expression for the transient correlator 
(|67p. Time-dependent projection operator manipulations 
combined with the theory of Volterra integral equations 
yield an exact equation of motion for the time evolu- 
tion of the transient correlator containing a generalized 
friction kernel - a memory function formed from the au- 
tocorrelation of fluctuating stresses. Mode-coupling type 
approximations to this kernel yield the nonlinear integro- 
differential equation p7l - [l9j 



(68) 



$q(i,t0) + rq(t,to)^q(i,to) 



+ J dt'm^{t,t',to)^q{t',to)j =0 

where the overdots denote partial differentiation with re- 
spect to the first time argument. Here the 'initial de- 
cay rate' obeys rq(t,to) = Doq'^{t,to) / Sq(t.to) with Dq 
the diffusion constant of an isolated particle. The for- 
mal manipulations presented in [itI [isj have revealed 
that imposing a time-dependent fiow results in a mem- 
ory kernel which depends upon three time arguments. 
The presence of a third time argument, which would have 
been difficult to guess on the basis of quiescent MOT in- 
tuitition, turns out to have important consequences for 
certain rapidly varying fiows (e.g. step strain [l3]) and 
is essential to obtain physically sensible results in such 
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cases. The memory kernel nin (t,t\ to) entering (|68)) is 
given by the factorized expression [1^ [l9| 



X yqkp(t', to) Kjkp(t, ^o)<i>k(t^t„)(^, t')^p(t'.to) (t, t') 

where p = q — k, and the vertex function obeys 



(69) 



Kjkp(i,to) = q(i,to) ■ (k(t,to)cfc(t,t^) + p(i,to)cp(t,to) 



(70) 



with Ornstcin-Zcrnicke direct correlation function ~ 
1 — 1/5*4; (sec Ea. (p9|) '). In the linear regime Eqs. (|66)) and 

reduce to the standard quiescent MCT forms 
and (dD), respectively. 

An important feature of Eos. ([70)1 is that they of- 
fer a closed constitutive equation requiring only the static 
structure factor and velocity gradient tensor K,{t) as in- 
put to calculating the stress tensor. The equilibrium 5*^ 
is determined by the interaction potential and thermo- 
dynamic statepoint and, as in quiescent MCT, serves as 
proxy for the pair potential (an interpretation arising 
from field-theoretical approaches to MCT j213l |). The 
role of Sq within the ITT-MCT should be contrasted 



with that within the Miyazaki-Reichman theory |201l - 
|203| |. discussed in subsection l7.11 where it enters as an ap- 
proximation to the flow-distorted structure factor S'(k). 

In subsection 12.31 we introduced the principle of mate- 
rial objectivity; an approximate symmetry requiring that 
a valid constitutive relation be rotationally invariant. 
While verification of rotational invariance is straightfor- 
ward for the phenomenlogical Lodge equation introduced 
in subsection l2.11 proof becomes more demanding for the 
microscopic constitutive theory given by Eas. (|66p - (|70p . 
Nevertheless, substitution of Eqs. ([TT|) and ([T^ into (|66p - 
([701) yields the resuU thus verifying that the ITT- 

MCT constitutive equation is indeed material objective 
as desired [l^ . Material objectivity is an important con- 
sistency check for constitutive theories based on over- 
damped Smoluchowski dynamics, for which it represents 
an exact symmetry constraint. 

Possibly the most exciting feature of the ITT-MCT 
constitutive equation ((66)) -([70 l) is that it incorporates a 
mechanism for describing the slow structural relaxation 
leading to dynamical arrest. The predicted rheological 
response thus goes from that of a viscous fluid to that 
of an amorphous solid, characterised by an elastic con- 
stant, upon variation of the thermodynamic control pa- 
rameters. The ability of the theory to unify the descrip- 
tion of fluid and glassy states stems from the underlying 
mode-coupling-type approximations which are tailored 
to capture the cooperative particle motion in dense col- 
loidal dispersions, ultimately leading to particle caging 
and arrest. Mathematically, this scenario arises from 
a bifurcation in the solution of the nonlinear integro- 
differential equation (|68|) at suffiently high volume frac- 
tion/attraction strength associated with a diverging re- 
laxation time of the transient density correlator. Glass 



formation within the MCT-ITT approach is a purely dy- 
namical phenomenon, as the equilibrium Sq used as in- 
put varies smoothly accross the transition. The fluid- 
solid transition contained within the ITT-MCT equations 
can be better appreciated by considering the small strain 
limit of Eq. (j66|) which yields the linear response result 

m 



\t)= dt' 



{(k.^(0-k)kk}^^^*^-(^-*') 



167r3 



kSk 



(71) 



where $fe(t) is the correlator from quiescent MCT jlOS 
In the glass the correlator does not relax to zero for long 
times and a partial integration of (|7ip followed by taking 
the limit of small strain leads to the result 



(T{t) = 2G{t^oo)e{t), 



(72) 



where e{t) is the infinitessimal strain tensor and G{t — > 
oo) is an elastic modulus obtained from Eq. p^ (also 
known as Lame's second coefficient, the first being zero 
here due to incompressibility) . Eq. ((7^ is essentially 
Hookes law, describing the small strain response of a 
glassy solid. Going beyond linear response, Eq. (j66|) incor- 
porates the fiuidizing effect of flow and thus makes pos- 
sible investigation of a large number of time-dependent 
rheological situations in which externally applied flow 
fields compete with glass formation and slow structural 
relaxation (see subsection 17. 6p . 



7.5. Distorted structure factor 

The microscopic ITT-MCT constitutive equation dis- 
cussed above enables comparisons to be made with tra- 
ditional continuum rheological modelling (section [2|) for 
which the macroscopic stress tensor is the fundamental 
quantity of interest. However, the formal ITT result ([55)1 
also enables calculation of the distorted structure fac- 
tor, S{'k,t) = I + nj dr{g{r,t) — l)e'^''^, which makes 
possible a comparison with the microstructure obtained 
from approaches based on the pair Smoluchowski equa- 
tion (section I5TT|) . In particular, setting / = Ap^^p^ = 
PkPk ^ (PtPk) "1 Eq. ([58)) yields the formal result 

S{k,t) = {plp^)+ f dt'{K{t'):aJj''''''^'^ApM. 



(73) 

Mode-coupling projection operator steps analagous to 
those leading to (|66p yield the ITT-MCT expression for 
the distorted structure factor 



dt' 



*k(^t')(t,t') (74) 



where the transient density correlator is given by solu- 
tion of Eas. ([55)) - (|70p for given Sk and K{t) and where an 
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isotropic term has been suppressed. Eq. ((74)) has the ap- 
peahng interpretation that flow-induced microstructural 
changes are built-up by integration of the affinely shifted 
equilibrium structure factor over the entire flow history, 
weighted by the transient density correlator describing 
the fading memory of the system. The temporally non- 
local character of Eg. ([7^ is in striking contrast to the 
local approximations based on the pair Smoluchowski 
equation. The former consists of a history integral over 
a memory function which is itself determined by solution 
of a nonlocal integro-difFerential equation whereas 
the latter are purely Markovian approximations. 

Within the pair Smoluchowski approach, in the ab- 
sence of HI the stress tensor is exactly related to the dis- 
torted pair correlation function by Eq. ([29| . It is therefore 
of interest to inquire whether a similar connection holds 
within the approximate ITT-MCT approach. It is a rel- 
atively simple exercise to show that, subject to a certain 
contraint to be discussed below, Eas. (p5)) and ([71)) are 
connected by the relation [l^ 



-ni - nknT 



dk kk 



16n 



3— c',<55kW, (75) 



where 6Su.it] k) = Sk{t;K) — Sk and 11 is the equilib- 
rium osmotic pressure. For shear flow <Jxy{t) from 
coincides with a result of Fredrickson and Larson [21 
for sheared copolymers, reflecting the Gaussian statis- 
tics underlying both the field-theory approach of [214 
and the ITT-MCT factorization approximations. Equa- 
tion ([75p thus connects stresses to microstructural distor- 
tions, which build up over time via the affine stretching 
of density fluctuations competing with structural rear- 
rangements encoded in $k(i,i')- 

For the off-diagonal stress tensor elements Eq. ([75|) con- 
nects Eas. ([55|) and ([71[) directly. For the diagonal ele- 
ments contibuting to the osmotic pressure Eq. ([75[) is also 
valid, providing that the following approximate 'sum- 
rule' is obeyed 



ad) Gtt 



dh\k 



d\x\6 



(76) 



Although it is not at all obvious that the above relation 
should hold, numerical calculations for hard spheres (us- 
ing e.g. the Percus-Yevick approximation for Sk |16| ) 
show that it represents a rather good approximation. 
It should be emphasized that the application of projec- 
tion operator approximations to (|65l) and (|73p yield ap- 
proximate expressions for the stress and structure factor, 
respectively, which are not neccessarily self-consistent, 
in the sense that integration of the approximate <S'(k, t) 
leads to the approximate cr{t). The fact that the ITT- 
MCT S'(k, t) is almost consistent with the direct ITT- 
MCT approximation to the stress is a testament to the 
underlying robustness of the method. 

Although the caging mechanism is expected to be most 
important for statepoints close to the glass transition, 
ITT-MCT calculations of S'(k) for hard-spheres under 



weak shear flow jl93[ | suggest that the (truncated) diver- 
gence of the structural relaxation timescale at the point 
of dynamical arrest (be it at the idealized glass transi- 
tion point, a s pre dicted by MCT, or some higher vol- 
ume fraction [iQ^) remains relevant for volume fractions 
well removed from the singularity and has a range of 
influence which extends back to dense equilibrium fluid 
states. More recently the MCT-ITT S'(k) has been eval- 
uated numerically for two-dimensional hard-discs under 
shear at finite values of Pe (79j . without invoking any 
additional isotropic approximations (see subsection 17. 6p . 
These calculations show only qualitative agreement with 
Brownian dynamics simulation results and overestimate 
the magnitude of the distortion by around an order of 
magnitude, a failing which is attributed to the fact that 
ITT-MCT apparently underestimates the speeding up of 
structural relaxations induced by the shear flow. This 
is to be contrasted with pair Smoluchowski-based ap- 
proaches (e.g. [l75, JL90J) which underestimate the mag- 
nitude of the low-shear distortion for dense fluid states 
{4> ~ 0.5) by around an order of magnitude. It there- 
fore appears that neither the pair Smoluchowski nor the 
ITT-MCT approach can account adequately for the shear 
induced distortion of the microstructure. 



7.6. Applications 

Explicit numerical solution of the ITT-MCT constitu- 
tive equation (Eas. ([^ - ([7n[) ) has been performed for a 
one-component system of hard spheres under steady flows 
of various geometry [isj . However, the computational re- 
sources required to solve the anisotropic Eqs. ([68| - ([70[) in 
three-dimensions over many decades of time are daunting 
and the efficient numerical algorithms required to reduce 
the computational load are still under development. Nev- 
ertheless, it is hoped that much of the essential physics 
may be captured by solving a simplified set of equations 
in which the advected wavevectors are approximated by 
an isotropic average k(i,i') — ?> kis{t,t'), where 

1 



kUt,t') 



in 



dnk^{t,t'). 



(77) 



This technical approximation has bee n su ccessfully ap- 
plied to the case of simple shear [sl, |203[ 1 and enables 
the angular integrals entering Eqs. ([68[) -([70 |) to be per- 
formed analytically. For two-dimensional systems algo- 
rithms have been developed which enable accurate nu- 
merical solution of the ITT-MC T eq uations without ad- 
ditional isotropic approximation [79[. 

In Fig(T3| fiow curves for hard-spheres resulting from 
solution of Eqs. ([55[) -([7(I [) subject to the isotropic approx- 
imation ([77|) are shown for both steady shear and steady 
planar extensional fiow. For these two choices of fiow the 
defining velocity gradient tensors are given by 



(78) 
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FIG. 14: Flow curves from the microscopic ITT-MCT ap- 
proach of [3 for hard-spheres at three different volume frac- 
tions close to the glass transition (j>c. Full lines show the 
steady shear stress, a = (Txy in units of kBT/{2Rf under 
shear flow as a function of Pe. Broken lines show the stress 
difference a = — Oyy (related to the extensional viscos- 
ity) under planar extensional flow. Each curve is labelled 
according to the distance in volume fraction from the glass 
transition, = (p — <j)c- (a) and (b) are fluid states with 
A0 = — 10~* and — 10"^, respectively. State (c) is in the glass, 
A(f> =10^*, and exhibits a dynamic yield stress for both flow 
geometries. The inset shows a possible realization of planar 
extensional flow. The lower two panels show the Trouton ra- 
tio (Txx — Cyy/cTxy as a function of Pe for both uniaxial and 
planar extensional flow. 



Specifically, FigfHl shows a^y under shear flow and Act = 
Cxx — CTyy under planar extension for hard-spheres as 
a function of Pe, for various volume fractions around 
the glass transition. The equilibrium structure factors 
used as input for these calculations were provided by the 
monodisperse Percus-Yevick theory [l^ . For extensional 
flows it is natural to plot the stress difference Act as a 
function of Pe, as this is simply related to the exten- 
sional viscosity rjc ~ Aa/j. Flow curves below the glass 
transition show a regime of linear response, character- 
ized by a constant viscosity, followed by shear thinning 
as Pe is increased. On approaching (f>g from below the 
linear regime moves to lower values of Pe and disappears 
entirely on crossing the glass transition. The resulting 
plateau in the flow curves identifies a dynamical yield 
stress (see subsection l5.4[) for both of the considered flow 
geometries. 

In the lower panels of Fig[T4] the Trouton ratio (ctxx — 
^yy)/'''xy is shown as a function of Pe for both pla- 
nar and uniaxial extensional flow fisj . For fluid states 
in the linear regime Trouton's rules assert that the ratio 



of extensional to shear viscosity ?7o/?7s takes the values 
4 and 3 for planar and uniaxial extension, respectively. 
These characteristic ratios arise from purely geometrical 
considerations and emerge naturally from the ITT-MCT 
approach as the Pe — >■ limiting values of the Trouton 
ratio. For glassy states (curves labelled (c) in Fig fH)) 
the structural relaxation time diverges and the linear re- 
sponse regime vanishes. As a consequence, the classical 
Trouton ratios arc not recovered for glassy states in the 
limit Pe — >■ and non-trivial values dictated by the dy- 
namical yield stress may be identifled. The results for ex- 
tensio nal fl ow presented in [l8| alongside those for simple 
shear |l40l | thus provide the first steps towards the predic- 
tion from first-prin ciple s theory of a dynamic yield stress 
surface for glasses [218j . Calculation of the yield surface 
from a simplified schematic version of Eqs. (|55)) - ((7n)) will 
be discussed in subsection 17.81 below. 

Going beyond steady flow, in [65j experiments on 
PMMA colloidal dispersions, molecular dynamics simu- 
lation and the ITT-MCT approach of were combined 
to study the evolution of stresses during start-up shear 
flow for high volume fraction fluids close to glassy ar- 
rest. The sudden onset of a steady shear flow leads to 
the build-up of stresses in the systems as a function of 
the accumulated strain 7 = jt. For small values of 7 
the response is elastic (as described by Eg. ([7^ ). whereas 
for large strains system enters steady state viscous flow 
with a stress independent of 7. In between these two lim- 
its, typically at strains around 10%, the stress exhibits 
a maximum as th e loc al microstructure is broken up by 
the external flow |l45l |. Although a stress 'overshoot' in 
response to start-up shear flow is a rather generic feature 
of the rhcology of complex fluids, its microscopic ori gins 
remain poorly understood. A central novel aspect of |65| 
was thus to connect the stress overshoot to anomolous 
behaviour in the mean-squared-displacement ('superdif- 
fusion') identified in both simulation and confocal mi- 
croscopy experiment. From the exact Eq. (j49p it is clear 
that the only way in which the shear stress can exhibit 
a maximum is if the modulus under steady shear, Gssit), 
becomes negative at long times. FigjS] shows that Gss{t) 
from the ITT-MCT approach (employing approximation 
([77)) ) indeed predicts negative values at long times. The 
inset to FiglS] shows the same data as a function of strain 
and demonstrates that the negative region of Gss{t) oc- 
curs at around 10% strain, consistent with the position 
of the stress overshoot. 



7.7. Schematic model 

The microscopic ITT-MCT constitutive equation out- 
lined in subsection 17.41 provides a route to first-principles 
prediction of the rheological behaviour of arrested col- 
loidal states. However, the anisotropic, wavevector de- 
pendent expressions are rather intractable for three di- 
mensional flows, hindering both their practical use and 
interpretation. In order to facilitate numerical calcula- 
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tions for flows of interest a simplified 'schematic' version 
of the tensorial microscopic theory has very recently been 
proposed [l^. Such schematic models have proved in- 
valuable in the analysis of mode-coupling theories and 
provide a simpler set of equations which retain the es- 
sent ial m athematical structure of the microscopic theory 

[Uliol. 

Applying the isotropic approximation (|77p to the mi- 
croscopic ITT-MCT expression for the stress ((66)) enables 
the angular integrals to be performed explicitly, leading 
to the simplified form 



ait) 



dt' 



d_ 

at' 



-B{t^t') 



G{t,t'), 



(79) 



where B is the Finger tensor and an explicit expression 
for G{t,t') may be found in [l^. By disregarding all 
wavevector dependence the modulus can be expressed in 
terms of a single-mode transient density correlator 



(80) 



where Va- = G{t, t) is a parameter measuring the strength 
of stress fluctuations (typically taking values of the order 
lOOksT / for hard-sphere-like colloids). A schematic 
equation-of-motion for ^(t,t') may be obtained by ne- 
glecting the wavevector dependence of the microscopic 
expression ([68)) . leading to 



(81) 

where the single decay rate F simply sets the time scale 
and may thus be set to unity. Experience with the con- 
struction of schematic MCT models both in the quiescent 
[l03l | and steady shear cases [S^ combined with analysis 
of the way in which strain enters the microscopic mem- 
ory (|69p lead to the following schematic ansatz for the 
memory function 

m{t,t',to) = h{t,to)h{t,t') [i^i^{t,t') + iy2<i>^{t,t')].{82) 

The parameters vi and i>2 represent in an unspecified 
way the role of Sq in the microscopic theory and, fol- 
lowing standard MCT practice, are given by V2 = 2 and 
vi = 2(V2 — 1) -t- e/{\/2 — 1). The separation parame- 
ter e encodes the distance from the glass transition, with 
negative values corresponding to fluid states and posi- 
tive values corresponding to glass states. Finally, the 
/i— function is given by 



72 + v[h {t, to) - 3) + (1 - i^){l2 (t, to) - 3) ■ 



(83) 



where 7c sets the strain scale (typically jc ~ 10%) and 
V is a mixing parameter (0 < < 1). The invariants of 
the Finger tensor, Ii = TrB and I2 = TrB"-*- incor- 
porate the fluidizing influence of flow into the memory 






FIG. 15: The state of stress of a material under applied 
force can be represented by a point in the three dimensional 
space of principal stresses. The cylinder shown here is the 
surface defined by the von Mises criterion (|86|l . Stress states 
lying within the cylinder are deformed by the applied stress, 
but do not yield, whereas states outside the cylinder exhibit 
plastic flow. Due to invariance along the 'hydrostatic' axis 
si — S2 = S3 (arising from incompressibility) the yield surface 
is usually viewed in the 'deviatoric' plane perpendicular to 
this (the lower left representation in this figure). 



function. Requiring that flow enter via the Finger tensor 
alone guarantees that the resulting schematic theory is 
material objective (see subsection 12. 3p , consistent with 
the fully microscopic theory. 

The schematic expression for the stress ((79)) is closely 
related to the Lodge equation of continuum mechanics 
(see subsection 12. ip . Integrating Ea. (|7^ by parts yields 



cT{t)^ 1* dt'B{t,t') (^-§p'^Ht,t') 



(84) 



Comparison of this expression with Eq.([5]) shows that 
the present schematic theory goes considerably beyond 
the standard Lodge equation by incorporating memory 
which is both nonexponential and a function of two time 
arguments, reflecting the loss of time-translational invari- 
ance under time-dependent flow. Replacing the correla- 
tor with a simple exponential trivially recovers the Lodge 
equation ([5]). 

The tensorial schematic theory given by Eqs. ((79)) - ((83)) 
has been applied to predict flow curves under shear and 
extensional flow, normal stresses under shear and the 
transient stress response to step strain. In all cases tested 
so far the predictions of the schematic model are in good 
qualitative agreement with those available from the mi- 
croscopic theory using the isot ropic approximation (|77p 
in three-dimensions [13, HI, Il40j | and exact numerical 
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FIG. 16: Left Panel: The dynamical yield surface from ITT- 
MCT [li] for a glass with e 3(0 - 0gt) = 10"^, where (jigt 
is the volume fraction at the hard-sphere glass transition, in 
the space of principal stresses (si, S2, S3) as viewed along the 
hydrostatic axis si — S2 = S3 (stress in units of / {2RY). 
The red points correspond to planar extensional flow and the 
blue points to uniaxial extensional flow. Right Panel: A closer 
view reveals that the surface is not a perfectly circular cylinder 
(as predicted by the von Mises criterion for static yielding) 
and that maximal deviation from circularity occurs at points 
of pure uniaxial extension. These deviations are connected 
with the existence of finite normal stress differences. 

solution in two-dimensions [t^. The schematic theory 
predicts a positive value for the first normal stress dif- 
ference Ni = (Txx — fyy undcr shear flow in accord with 
microscopic ITT-MCT calculations. On the other hand, 
the second normal stress difference N2 = (Jyy^<J2z is from 
Ea.(|79p identically zero, in disagreement with both an- 
alyti cal l ow Pe analysis of the pair Smoluchowski equa- 
tion |l6l| and colloidal experiments j215{ (finite negative 
values are found). The disappearance of N2 from the 
schematic theory can be traced back to the isotropic ap- 
proximation leading to Eg. ((7^ which effectively kills off 
this feature of the fully microscopic theory. 

7.8. Yield stress surface 

A striking feature of the theory developed in [l^ is that 
it permits direct calculation of a dynamic yield stress sur- 
face for glasses. Related static yield surfaces have been 
empirically postulated and employed for over a century 
in the engineering community to study the yielding of 
amorphous solids (see also subsection l5.4p . The two clas- 
sical criteria for deter mining the onse t of plastic yield 
are due to Tresca |216j and von Mises j217j . The Tresca 
criterion asserts that a material will yield when the max- 
imum shear stress due to the deformation exceeds a crit- 
ical value. Recalling that an external force imposed on 
a material can be represented as a stress tensor which 
can be diagonalized to obtain values for the three prin- 
cipal stresses si,S2 and S3, t he T resca criterion can be 
compactly stated in the form |218l | 

Max( I si - S2I, |S2 - S3U si - sal ) = ^Piul^ , (85) 

where ^ is the shear stress at yield under a simple shear 
deformation. According to Ea. ([551) knowledge of is 



thus sufficient to determine the mechanical stability of a 
material under an arbitrary applied force. Alternatively, 
the von Mises critierion requires that th e dist ortion strain 
energy exceeds a critical value at yield j218l | 

\ ((51 - s^f + (S2 - s^f + (si - S3)2) = {al,f . (86) 

Both the Tresca and von Mises criteria have proven to 
be in reasonable qualitative agre ement with yield experi- 
ments on crystalline metals |219{ . Two main assumptions 
undcrly Eqs. ()85|) and ()86|) : (i) the microscopic rearrange- 
ments leading to plastic deformation do not lead to sig- 
nificant dilation of the material, (ii) that residual stresses 
arising from the deformation history of the sample do not 
influence the yielding (i.e. there is no Bauschinger effect). 

It is useful to interpret equation ()86p geometrically in 
the space of principal stresses, where it describes a sur- 
face separating elastically deformed states from states of 
plastic flow. Equation ([SS)) defines a circular cylinder 
with axis along the line si = S2 = S3 and radius v^cf^,. 
The symmetry about this 'hydrostatic' axis is a geomet- 
rical reflection of the fact that the yield condition ([55)1 
is independent of hydrostatic pressure. The plane which 
passes through the origin and which lies perpendicular 
to the cylinder axis is the so-called deviatoric plane. All 
yield stress surfaces which are independent of hydrostatic 
pressure may be projected without loss of information 
onto the deviatoric plane. In the case of the von Mises 
and Tresca criteria this generates a circle (see Fig fTS]) and 
a hexagon, respectively. 

The results presented in [l9| reveal intriguing connec- 
tions between the static yielding discussed above and dy- 
namic yield, as determined from the Pe — >■ values of 




FIG. 17: The yield stress surface in the t^b (rotational shear 
stress) Trz (squeeze shear stress) plane, where the stresses have 
been scaled by the yield stress Tc in simple shear. The blue 
and red squares show results for two different emulsions with 
differing values of Tc (28Pa and 52Pa). The green squares 
are data taken using a carbopol gel (tc — 70Pa). In this rep- 
resentation the von Mises criterion (|86[l becomes a circle and 
is indicated by the solid line. Reprinted by permission from 
MacMillan Publishers Ltd: Nature Materials 9 115, copyright 
(2010). 
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the flow curves. Within the ITT-MCT approach, for any 
given steady flow field (e.g. shear, uniaxial extension) 
there exists for glassy states a finite stress tensor in the 
limit of vanishing flow rate. This stress tensor at yield 
may be diagonalized to obtain three eigenvalues and plot- 
ted as a point in the cartesian space of principal stresses. 
By considering all possible nondcgcncratc flows (made 
possible by a suitable parameterization of the velocity 
gradient tensor [l^ ) a closed locus of points may be con- 
structed in the space: The dynamic yield stress surface. 
In [l^ this procedure was realized using the schematic 
model Eqs.([7gi)-(P1). 

In Fig llGI we show the deviatoric projection of the dy- 
namical yield stress surface of a colloidal glass, as pre- 
dicted by the schematic model. To a first approximation 
the numerically calculated dynamical yield surface from 
the theory agrees well with the empirical von Miscs crite- 
rion for static yield However, closer inspection (see 
FiglTSl right panel) reveals that deviations at around the 
percent level occur, with the maximal deviation located 
at points of pure uniaxial extension. A careful analysis 
of the schematic equations reveals that this fine struc- 
ture of the surface can be attributed to the existence 
of a finite first normal stress difference [l^. Expanding 
the schematic model result to first order in iVf , the first 
normal stress difference at yield, provides an explicit ex- 
pression for the schematic model yield surface 



1 



12^ 1^ + (3 + ^2)3/2''!^-' ^^'^ 



where < A < 1 parameterizes the geometry of the im- 
posed flow. Eq. (j87|) thus describes a noncircular cylinder 
with a radius which varies according to the value of the 
parameter A, each value of which corresponds to a given 
azimuthal angle about the hydrostatic axis. Higher order 
terms in the expansion ()87|) exist (and can be explicitly 
calculated) but remain numerically negligable due to the 
smallness of Ni/a^^. 

Very recent experiments on yielding soft materials have 
attempted to determine the shape of th e yie ld surface 
[220l |. The novel rheometer employed in |22Cll | applies a 
combined squeeze and rotational shear flow to a material 
sample loaded between two parallel plates. The claim is 
that by independently varying the rotation and squeeze 
rate it may become possible to explore the entire family 
of flows in a way analogous to the mathematical param- 
eterization of the velocity gradient tensor employed in 
[l9| to calculate the schematic model yield surface. How- 
ever, it remains uncertain whether the superposition of 
two shear flows (radial Poiscuille and tangential homoge- 
neous shear) is really sufficient to map the yield surface. 
Fig IT71 shows experimental data taken using three differ- 
ent yieldin g ma terials: a Carbopol gel and two different 
emulsions [2201 ■ 

The yield data are shown in the rotational shear stress 
(Ti-e), squeeze shear stress (xrz) plane. While the data 



shown in Fig ) 171 are not inconsistent with the von Mises 
criteria, it may well be that the chosen superposition flow 
actually constrains the surface to be spherical, regardless 
of the true form of the yield surface. It would cert ainly 
be remarkable if the soft materials considered in |220l | 
obey yield criteria developed for crystalline solids, de- 
spite the ver y dif ferent underlying microscopic plasticity 
mechanisms 
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Whether the microscopic theory (Eas. ([S5)) - (|7n)) ) in- 
deed predicts a similar dynamic yield surface, as ex- 
pected, and its relationship to static yielding in glassy 
materials remain important open problems. Neverthe- 
less, the results from the schematic model are promising 
and represent a considerable step towards a microscopic 
derivation of material specific yield surfaces from first 
principles. 



8. OUTLOOK 

In this review we have attempted to provide an 
overview of the rheological phenomenology presented by 
colloidal dispersions and to outline some of the leading 
theoretical approaches aiming to rationalize this. It is 
clear that much remains to be done and that cxisiting 
theories have met with only partial success in solving the 
complex many-body problem of driven, strongly inter- 
acting colloidal systems. We hope that both the presen- 
tation and choice of topics contained within the present 
work serve to emphasize the common ground between dif- 
ferent theories (continuum mechanics, pair Smoluchowski 
treatments and Green-Kubo based approaches), as well 
as to highlight where progress still needs to be made. 

One of the clear deficiencies of approximate closures 
of the pair Smoluchowski equation is their apparent in- 
ability to describe, even qualitatively, the slow structural 
relaxation time present in colloidal dispersions at finite 
volume fraction. This failing is inherent in the Marko- 
vian character of the approximate closures, which follows 
as a natural consequence of applying equilibrium statisit- 
ical mechanical relations to connect pair and triplet cor- 
relation functions. A clear challenge to future theo- 
ries which attempt to improve this situation is thus to 
tackle directly the intinsic difference between nonequlib- 
rium and equlibrium by confronting the irreducible term 
(see Ea l46)) containing the missing physics. Moreover, 
the existing numerical data in the literature is restricted, 
largely for technical reasons, to the low Pe limit. While 
this enables interesting analysis of the zero-shear viscos- 
ity and leading order microstructural distortion, there is 
an absence of data for the nonlinear rheology as a func- 
tion of Pe. 

Despite the success of the ITT-MCT approach in de- 
scribing the nonlinear response of states close to dynamic 
arrest, there remain gaps in the theoretical formulation 
which should be filled and many fundamental questions 
are yet to be addressed. A notable omission is that the 
present formulation of the theory docs not enable incor- 



36 



poration of hydrodynamic interactions. The influence of 
hydrodynamics on the rheology of densely packed glassy 
states is largely unexplored and their incorporation into 
the theory, even at a crude approximate level, would 
therefore be of considerable interest. Moreover, due to 
the recent nature of the ITT-MCT theory and the com- 
plexity of numerically solving the equations, many flows 
of rheological interest remain to be explored in detail. 
Of particular importance is to assess the predictions of 
the theory for nonsteady flows where interesting inter- 
action effects between yielding and the time-dependent 
strain field may be envisaged (e.g. large amplitude oscil- 



[1] R. Mezzenga, P. Schurtenberger, A. Burbidge and 
M. Michel, Nature Materials 4 729 (2005). 

[2] R. G. Larson, The structure and rheology of complex 
fluids (Oxford University Press, New York, 1998). 

[3] E. Brown and H.J. Jaeger, Phys.Rev.Lctt. 103 086001 
(2009). 

[4] A. Fall, N. Huang, F. Bertrand, G. Ovarlez and 
D. Bonn, Phys.Rev.Lett. 100 018301 (2008). 

[5] P. Coussot, Rheometry of pastes, suspensions and gran- 
ular materials (Wiley, 2005). 

[6] W.B. Russel, D.A. Saville and W.R. Schowalter, Col- 
loidal dispersions (Cambridge University Press, 1989). 

[7] W. van Megen and P. Pusey, Phys.Rev.A 43 5429 

(1991) . 

[8] W. van Megen and S. Underwood, Phys.Rev.Lett. 70 
2766 (1993). 

[9] W. van Megen and S. Underwood, Phys.Rev.E 49 4206 
(1994). 

[10] For example, in Refs.jg] and [y| dynamic light scatter- 
ing was employed to measure the coherent transient 
density correlator of spherical hard-sphere- like colloids. 
The data confirmed the factorization properties of the 
a and /3 relaxation, as predicted by the mode-coupling 
theory of the glass transition [ll|, [l^l. In [l^ confocal 
microscopy experiments on a mixture of PMMA col- 
loids and non-adsorbing polymer were used to confirm 
the capillary wave theory of the fluctuating interface 
between demixed fluid phases [T^ . 

[11] W. Gotze and L. Sjogren, Rep.Prog.Phys. 55 241 

(1992) . 

[12] W. Gotze, J.Phys.:Condens.Matter 11 Al (1999). 

[13] D.G.A.L. Aarts et al, science 304 847 (2004). 

[14] L. Mandelstam, Ann.Phys. 41 609 (1913). 

[15] J. K. G. Dhont, An introduction to the dynamics of col- 
loids (Amsterdam, Elsevier, 1996). 

[16] J. -P. Hansen and I.R. McDonald. Theory of Simple Liq- 
uids (Academic, London, 1986). 

[17] J.M. Brader, Th. Voigtmann, M.E. Gates and M. Fuchs, 
Phys.Rev.Lett. 98, 058301 (2007). 

[18] J.M. Brader, M.E. Gates and M. Fuchs, Phys.Rev.Lett. 
101 138301 (2008). 

[19] J.M. Brader et al, Proc. Natl. Acad. Sci. U.S.A. 106 
15186 (2009). 

[20] M. Fuchs, in High solid dispersions, ed. M. Cloitre 



latory shear, where higher harmonics will contribute to 
the stress response). 



Acknowledgments 

I would particularly like to thank M. Fuchs, M.E. 
Gates, Th. Voigtmann and M. Kriiger for many stim- 
ulating discussions. Support was provided by the SFB 
TR6 and the Swiss National Science Foundation. 



(Springer, Berlin, 2009). 

[21] M. Fuchs and M.E. Gates, J.Rheol. 53 957 (2009). 

[22] M. Doi and S. F. Edwards, The Theory of Polymer Dy- 
namics (Oxford University Press, 1989). 

[23] R.G. Larson, Constitutive equations for polymer melts 
and solutions Butterworth, Boston, 1988). 

[24] T.G.B. Macleish, Adv.Phys. 51 1379 (2002). 

[25] M.E. Gates and S.M. Fielding, Adv.Phys. 55 799 (2006). 

[26] P. Coussot, Soft Matter 3 528 (2007). 

[27] J. Vermant and M.J. Solomon, J. Phys.:Gondens. Matter 
17 R187 (2005). 

[28] R.A. Lionberger and W.B. Russel, Adv.Chem.Phys. 
Ill 399 (2000). 

[29] J.F. Morris, Rheol.Acta 48 827 (2009). 

[30] R.B. Bird, R.G. Armstrong and O. Hassager, Dynamics 
of polymeric liquids (Wiley, 1987). 

[31] G.K. Batchelor, An Introduction to fluid dynamics 
(Cambridge University press, 1967). 

[32] G. Truesdell and W. Noll, The nonlinear field theories 
of mechanics (Springer, Berlin, 1965). 

[33] A useful discussion of material objectivity may be found 
in [s^j. In Addition, we refer the reader to Physica A 
118 43 (1983), which documents the insightful com- 
ments of P.G. de Gennes regarding this issue. 

[34] G. Ryskin, Phys.Rev.A 32 1239 (1985). 

[35] P. SoUich, F. Lequeux, P. Hebraud and M.E. Gates, 
Phys.Rev.Lett. 78 2020 (1997). 

[36] P. SoUich, Phys.Rev.E 58 738 (1998). 

[37] S.M. Fielding, P. SoUich and M.E. Gates, J.Rheol. 44 
323 (2000). 

[38] J. Crocker and D.G. Crier, J. Colloid Interface Sci. 179 
298 (1996). 

[39] V. Prasad, D. Semwogerere and E.R. Weeks, 
J. Phys.Cond. Matter 19 113102 (2007). 

[40] R. Besseling, L. Isa, E.R. Weeks and W.C.K. Poon, Ad- 
vances in Col. Int. Science 146 1 (2009). 

[41] F. Varnik, J.Chem.Phys. 125 164514 (2006). 

[42] A. Sierou and J.F. Brady, J. Fluid Mech. 448 115 

(2001) . 

[43] B. Doliwa and A. Heuer, Phys.Rev.E 61 6898 (2000). 
[44] L. Berthier, J.-L. Barrat and J. Kurchan, Phys.Rev.E 

61 5464 (2000). 
[45] L. Berthier and J.-L. Barrat, J.Chem.Phys. 116 6228 

(2002) . 

[46] J. -P. Bouchaud, F. Gugliandolo, J. Kurchan and 
M. Mezard, in Spin glasses and random fields, ed. 



37 



A.P. Young (World Scientific, Singapore, 1998). 
[47] H. Riskin, The Fokker-Planck equation (Springer, [71] 
Berlin, 1989). 

[48] S. Kim and S.J. Karilla, Microhydrodynamics, princi- [72 

pies and selected applications (Butterworth-Heinemann, [73 

Boston, 1991). [74 
[49] H. Jorquera and J.S. Dahler, J.Chem.Phys. 101 1392 

(1994). [75 
[50] J.F. Brady and G. Bossis, Ann.Rev.Fluid Mech. 20 111 

(1988). [76 
[51] A.J. Banchio and J.F. Brady, J.Chem.Phys. 118 10323 

(2003). [77 
[52] J.F. Brady, Chemical Engineering Science 56 2921 

(2001). [78 
[53] S. Sami, Stokesian dynamics simulation of extensional [79 

flow of Brownian suspensions, 

(Masters Thesis, California Institute of Technology, [80 
1997). 

[54] J. Mathews and R.L. Walker, Mathematical methods of [81 

physics (Addison- Wesley, 1970). 
[55] R.T. Foister and T.G.M. van de Ven, J.Fluid.Mech. 96 [82 

105 (1980). 

[56] G.I. Taylor, Proc.Roy.Soc.A 219 186 (1953). [83 
[57] J.F. Brady and J.F. Morris, J. Fluid Mech. 348 103 [84 

(1997). [85 
[58] M. Fuchs and M.E. Gates, Faraday Discussions 123, 267 [86 

(2003). 

[59] J.W. Bender and N.J. Wagner, J. Rheol. 40, 899 (1996). [87 

[60] Interestingly, the force acting on an isolated hard-sphere [i 
of radius R in a Newtonian solvent under (Stokesian) [89 
shear flow is identical to that acting on a point particle 
(i? — ^ 0) moving with the affine flow; i.e. the fact that [90 
the solvent flows around the particle does not influence [91 
the resulting force. However, setting v(r, t) = K,{t) ■ r is 
an approximation for finite colloidal volume fraction. [92 

[61] In [63 ] a confocal microscopy study of polymethyl- 
methacrylate (PMMA) colloids suspended in a mixture [93 
of decalin and tetralin revealed inhomogeneous plug flow 
for glassy states. In contrast, the rheological experi- [94 
ments performed in [6^ and [M ] using thermosensitive 
PNIPAM core-shell particles are fully consistent with [95 
homogeneous flow, both above and below the glass tran- 
sition. [96 

[62] P. Ballesta, R. Besseling, L. Isa, G. Petekidis and 

W.C.K. Poon, Phys.Rev.Lett. 101 258301 (2008). [97 

[63] J.J. Crassous et al, J.Chem.Phys. 128 204902 (2008). 

[64] M. Siebenbiirger, M. Fuchs, H. Winter and M. Ballauff, 

J.Rheol. 53 707 (2009). [98 

[65] J. Zausch et al, J. Phys.:Condens. Matter 20 404210 [99 
(2008). 

[66] M. Rauscher, A. Dominguez, M. Kriiger and F. Penna, 

J.Chem.Phys. 127 244906 (2007) 
[67] B.J. Alder and T.E. Wainwright, J.Chem.Phys. 27 1208 

(1957), W.G. Hoover and F.H. Ree, J.Chem.Phys. 49 

3609 (1968). 

[68] U. Gasser, J.Phys.Condens.Matter 21 203101 (2009). 

[69] In order to observe glassy behaviour the timescale of 

crystal nucleation Tfr has to be extended beyond the ac- [100] 
cessible experimental or simulation measurement time. 
Fortunately, rfr is a sensitive function of the degree of [101] 
polydispersity and for hard-sphere-like systems a poly- 
dispersity ~ 5% is usually sufficient to effectively su- [102] 
press crystallization. [103] 

[70] Th. Voigtmann, A.M. Puertas and M. Fuchs, 



Phys.Rev.E 70 061506 (2004). 

M.E. Gates and S. Ramaswamy, Phys.Rev.Lett. 96 
135701 (2006). 

S-H. Chong, Phys.Rev.E 78 041501 (2008). 

L. Berthier et al. Science 310 5755 (2006). 

G. Biroli, J-P. Bouchaud, K. Miyazaki and D.R. Reich- 

mann, Phys.Rev.Lett. 97 195701 (2006). 

A. Duri, H. Bissig, V. Trappe and L. Cipelletti, 

Phys.Rev.E 72 051401 (2005). 

R. Cerbino and V. Trappe, Phys.Rev.Lett. 100 188102 
(2008). 

E. R. Weeks, J.C. Crocker, A.C. Levitt, A. Schofield and 
D.A. Weitz, Science 287 627 (2000). 

M. Bayer, et al. Phys.Rev.E 76 011508 (2007). 

O. Henrich, F. Weyfier, M.E. Gates and M. Fuchs, 

Phil.Trans.R.Soc.A 367 5033 (2009). 

W. Gotze and Th. Voigtmann, Phys.Rev.E 67 021507 

(2003). 

D. Hajnal, J.M. Brader and R. Schilling, Phys.Rev.E 
80 021503 (2009). 

V. Trappe and P. Sandkuhler, Current Opinion in Col- 
loid & Interface Science 8 494 (2004)). 
P.N. Segre et al. Phys.Rev.Lett. 86 6042 (2001). 
K.N.J. Pham et al. Science 296 104 (2002). 
K.N.J. Pham et al. Phys.Rev.E 69 011503 (2004). 
M. Dijkstra, J.M. Brader and R. Evans, 
J.Phys.:Condens.Matter 11 10079 (1999). 
J.M. Brader and R. Evans, Physica A 306 287 (2002). 
M. Laurati et al, J.Chem.Phys. 130 134907 (2009). 
C.J. Dibble, M. Kogan and M.J. Solomon, Phys.Rev.E 

74 041403 (2006). 

S. Mauley et al, Phys.Rev.Lett. 94 238302 (2005). 

G. Foffi, C. De Michele, F. Sciortino and P. Tartaglia, 
Phys.Rev.Lett. 94 078301 (2005). 

F. Sciortino, P. Tartaglia and E. Zaccarelli, 
J.Phys.Chem.B 109 21942 (2005). 

A.J. Archer and N.B. Wilding, Phys. Rev. E 76 0315011 
(2007). 

C.L. Klix, CP. Royall, and Hajime Tanaka, Phys. Rev. 
Lett. 104, 165702 (2010). 

E. Bianchi, J. Largo, P. Tartaglia, E. Zaccarelli and 

F. Sciortino, Phys.Rev.Lett. 97 168301 (2007). 

E. Zaccarelh, J. Phys. :Condens. Matter 19 323101 
(2007). 

For example, the yielding behaviour of attractive col- 
loids is considered in K.N. Pham, et al Europhys. Lett. 

75 624 (2006). 

H. M. Laun, Angew.Makromol.Chem. 123 335 (1984). 
In general, the initial value of the shear modulus is de- 
termined by near-field hydrodynamic lubrication forces. 
The commonly studied case of Brownian hard-spheres 
in the absence of hydrodynamic interactions is a patho- 
logical special case for which G{t 0) ~ t^a as a 
consequence of the discontinuous potential interaction. 
Nevertheless, the Brownian hard-sphere G{t) rapidly 
be come s physical for later times. We refer the reader 
to [lOOl ] for more details on this issue. 

R.A. Lionberger and W.B. Russel, J.Rheol. 38 1885 
(1994). 

G. Nagele and J. Bergenholtz, J.Chem.Phys. 108 9893 
(1998). 

Z. Cheng et al, Phys. Rev. E 65, 041405 (2002). 

W. Gotze, in Liquids, Freezing and the Glass Transi- 



38 



J.Rheol. 36 743 (1992). 
Bung and F. Schmidt, J.Rheol. 35 999 

.F. Brady, J.Chem.Phys. 91 1866 (1989). 
Wagner, Rheologica Acta, 49 897 



tion ed. J-P. Hansen, D. Levesque and J. Zinn- Justin 
(Amsterdam: North-Holland, 1991). 
[104] T. Phung and J.F. Brady, in Slow Dynamics in 
Condensed Matter, K. Kawasaki, M. Tokuyama and 
T. Kawakatsu, (AIP, Woodbury, NY, 1992) Vol.256 
p.391. 

[105] S.R. Rastogi and N.J. Wagner, Comp.Chem.Eng. 19 
693 (1995). 

[106] G. Bossis and J.F. Brady, J.Chem.Phys. 80 5141 (1984). 
[107] J.R. Melrose and R.C. Ball, J.Rheol. 48 937 (2004). 
[108] J.R. Melrose and R.C. Ball, J.Rheol. 48 961 (2004). 
[109] S.D. Kulkarni and J.F. Morris, J.Rheol. 53 417 (2009). 
[110] J. Mewis and N.J. Wagner, Adv. Colloid Interface Sci. 

147-148 214 (2009). 
[Ill] P. D'Haene, J. Mewis and G.G. Fuller, J.CoU.Int.Sci. 

156 350 (1993). 
[112] H.A. Barnes, J.Rheol. 33 329 (1989). 
[113] R.L. Hoffman, Trans. Soc.Rheol. 16 155 (1972). 
[114] W.H. Boersma, J. Laven and H.N. Stein, AIChE J. 36 

321 (1990). 
[115] H.M. Laun et 
[116] H.M. Laun, R. 

(1991). 
[117] G. Bossis and J. 
[118] D.P. Kalman and N.J. 

(2009). 

[119] R.L. Hoffman, J.Rheol. 42 111 (1998). 
[120] R.C. Egres and N.J. Wagner, J.Rheol. 49 719 (2005). 
[121] H. Laun, J.Non- Newtonian Fluid Mech. 54 87 (1994). 
[122] W.J. Frith, P. d'Heane, R. BuscaU and J. Mewis, 

J.Rheol. 40 531 (1996). 
[123] M.E. Gates, C.B. Holmes, M. Fuchs and O. Henrich, in 

Unifying concepts in granular media and glasses, edited 

by A. Coniglio et al. (Elsevier, Amsterdam, 2004). 
[124] C.B. Holmes, M.E. Gates, M. Fuchs and P. SoUich, 

J.Rheol. 49 237 (2005). 
[125] M.E. Gates, M.D. Haw and C.B. 

J. Phys.:Condens. Matter 17 S2517 (2005). 
[126] J. Delhommelle and J. Petravic, J.Chem.Phys. 

074707 (2005). 

[127] A. Mehta, G. Barker and J.M. Luck, Physics Today, 62 
(5) 40 (2009). 

[128] V. Gopalakrishnan and C.F. Zukoski, J.Rheol. 48 1321 
(2004). 

[129] A.S. Negi and CO. Osuji, Rheologica Acta, 48 871 
(2009). 

[130] W. Poon, Science 304 830 (2004). 

[131] D.J. Evans and G.P. Morriss, Phys.Rev.Lett. 56 2172 
(1986). 

[132] D.J. Evans et al, Phys.Rev.A 46 6731 (1992). 

[133] M.L. Matin, P.J. Daivis and B.D. Todd, J.Chem.Phys. 
113 9122 (2000). 

[134] P.J. Daivis, M.L. Matin and B.D. Todd, J. Non- 
Newtonian Fluid Mech. Ill 1 (2003). 

[135] D. Lootens, H. van Damme, Y. Hemar and P. Hebraud, 
Phys.Rev.Lett. 95 268302 (2005). 

[136] P. Hebraud, Rheologica Acta 48 845 (2009). 

[137] D.R. Foss and J.F. Brady, J. Rheol. 44 629 (2000). 

[138] J. Bergenholtz, J.F. Brady and M. Vicic, J.Fluid.Mech. 
456 239 (2002). 

[139] M. Zackrisson, A. Stradner, P. Schurtenberger and J. 
Bergenholtz, Phys.Rev.E 73 011408 (2006). 

[140] M. Fuchs and M.E. Gates, Phys.Rev.Lett. 89 248304 
(2002). 



[141] 
[142] 



Holmes, 



123 



[143] 



144 

145; 

146; 

147 

148; 

149' 
150^ 
151 
152 
153^ 
154' 
155; 

156; 

157 

158; 

159; 

160 
161 
162; 

163 
164 
165' 



166; 

167 
168; 

169' 
170' 
171 
172 
173 

174 
175' 
176 
177 



R. Besseling, E.R. Weeks, A.B. Schofield and 
W.C.K. Poon, Phys.Rev.Lett. 99 028301 (2007). 
We note that there exist numerous alternative defini- 
tions of the 'yield stress' in the literature. Two common 
choices are: (i) the maximum value of the stress over- 
shoot which occurs in response to the onset of steady 
shear flow, (ii) the point for which G' = G" as a func- 
tion of amplitude for fixed frequency in oscillatory ex- 
periments (strain sweeps). 

This expecation is supported by the simulations per- 
formed in [I45I ] which found that the maximum of 
the stress-strain curve during start-up shear flow, often 
identified as the static yield-stress, exhibits a logarith- 
mic dependence on the age of the system. 
G. Ovarlez, S. Rodts, X. Chateau and P. Coussot, Rhe- 
ologica Acta, 48 831 (2009). 

F. Varnik, L. Bocquet and J.-L. Barrat, J.Chem.Phys. 
120 2788 (2004). 

F. Varnik, L. Bocquet, J.-L. Barrat and L. Berthier, 
Phys.Rev.Lett. 90 095702 (2003). 

N.J. Wagner and W.B. Russel, Physica A 155 475 
(1989). 

W.B. Russel and A.P. Cast, J.Chem.Phys. 84 1815 
(1985). 

N.J. Wagner, J. Colloid Interface Sci. 161 169 (1994). 
S. Kim and R.T. Mifhin, Phys.Fluids 28 2033 (1985). 
M. Medina-Noyola, Phys.Rev.Lett. 60 2705 (1988). 
J.F. Brady, J.Chem.Phys. 99 567 (1993). 
P. Attard, J.Chem.Phys. 91 3072 (1989). 
J.M. Brader, J.Chem.Phys. 128 104503 (2008). 
J.G. Kirkwood, F.P. Buff and M.S. Green, 
J.Chem.Phys. 17 988 (1949). 

Y. Yurkovetsky and J.F. Morris, J.Chem.Phys. 124 
204908 (2006). 

A. Einstein, Ann. Phys. 324, 289 (1906), see also the 
erratum A. Einstein, Ann. Phys. 339, 591 (1911). 

G. K. Batchelor and J.T. Green, J. Fluid Mech. 56 375 
(1972). 

G.K. Batchelor and J.T. Green, J. Fluid Mech. 56 401 
(1972). 

G.K. Batchelor, J.Fluid.Mech. 83 97 (1977). 
J.F. Brady and M. Vicic, J.Rheol. 39 545 (1995). 

F. Gadala-Maria and A. Acrivos, J.Rheol. 24 799 
(1980). 

J.K.G. Dhont, J.Fluid.Mech. 204 421 (1989). 
R.A. Lionberger, J.Rheol. 42 843 (1998). 
S. Ramaswamy, in Theoretical challenges in the dynam- 
ics of colloidal fluids edited by T. McLeish (Kluwer, 
Dordrecht, 1997). 

J. Blawzdziewicz and G. Szamel, Phys.Rev.E 48 4632 
(1993). 

J.F. Brady, J.Fluid.Mech. 272 109 (1994). 

M. Born and M.S. Green, A General Kinetic Theory of 

Liquids (Cambridge University Press, 1949). 

J.G. Kirkwood, J.Chem.Phys. 3 300 (1935). 

S. Rice and J. Lekner, J.Chem.Phys. 42 3559 (1965). 

T. Ohtsuki, Physica A 108A 441 (1981). 

T. Ohtsuki, Physica A 110 A 606 (1982). 

P. Attard and G. Stell, Chem.Phys.Lett. 189 128 

(1992). 

B. Cichocki, Z.Phys.B 66 537 (1987). 

G. Szamel, J.Chem.Phys. 114 8708 (2001). 

L. Verlet and J.-J. Weis, Phys.Rev.A 5 939 (1972). 
G. Stell, in Equilibrium theory of classical fluids (H.L. 



39 



Frisch and J.L. Lebowitz eds.) Benjamin, New York 
(1964). 

[178] J.M. Brader, Int.J.Thermophys. 27 394 (2006). 
[179] L. Verlet, Physica 30 95 (1964). 

[180] D. Levesque and L. Verlet, Physics Letters 11 36 (1964). 

[181] L. Verlet, Physica 31 959 (1965). 

[182] L. Verlet, Physica 32 304 (1966). 

[183] L. Verlet and D. Levesque, Physica 254 254 (1967). 

[184] M.S. Wertheim, J.Math.Phys. 8 927 (1967). 

[185] R.J. Baxter, Annals of Physics 46 509 (1968). 

[186] G. Stell, Mol.Phys. 16 209 (1969). 

[187] J-L. Barrat, J-P. Hansen and G. Pastore, 

Phys.Rev.Lett. 58 2075 (1987). 
[188] P. Attard, J.Chem.Phys. 91 3072 (1989). 
[189] K. Scherwinski, Mol.Phys. 70 797 (1990). 
[190] R.A. Lionberger and W.B. Russel, J.Chem.Phys. 106 

402 (1997). 

[191] F.J. Rogers and D.A. Young, Phys.Rev.A 30 999 (1984). 
[192] R.A. Lionberger and W.B. Russel, J.Rheol. 399 41 
(1997). 

[193] O. Henrich, O. Pfeifroth and M. Fuchs, 
J.Phys.:Condens. Matter 19 205132 (2007). 

[194] P. Resibois and J.L. Lebowitz, J.Stat.Phys. 12 483 
(1975). 

[195] D. Ronis, Phys.Rev.A 29 1453 (1984). 

[196] J.F. Schwarzi and S. Hess, Phys.Rev.A 33 4277 (1986). 

[197] N.J. Wagner, J. Non-Newtonian Flui Mech. 96 177 

(2001) . 

[198] H.C. Ottinger, Beyond equilibrium thermodynamics 

(Wiley, Hoboken, 2005). 
[199] D.J. Evans and G. Morriss, Statistical mechanics of 

nonequilibrium liquids, (Cambridge University Press, 

2008). 

[200] A.V. Indrani and S. Ramaswamy, Phys.Rev.E 52 6492 
(1995). 

[201] K. Miyazaki and D.R. Reichman, Phys.Rev.E 66 050501 

(2002) . 

[202] K. Miyazaki, D.R. Reichman and R. Yamamoto, 

Phys.Rev.E 70 011501 (2004). 
[203] K. Miyazaki, H.M. Wyss, D.A. Weitz and D.R. Reich- 

mann, Europhys.Lett. 75 915 (2006). 



[204] V. Koblev and K.S. Schweizer, Phys.Rev.E 71 021401 
(2005). 

[205] E.J. Saltzman, G. Yatsenko and K.S. Schweizer, 
J. Phys.:Condens. Matter 20 244129 (2008). 

[206] G. Petekidis, D. Vlassopoulos and P.N. Pusey, 
J. Phys.:Condens. Matter 16 S3955 (2004). 

[207] M. Fuchs and M.E. Gates, J.Phys.:Gondens.Matter 15 
S401 (2003). 

[208] M. Fuchs and M. Ballauff, J.Ghem.Phys. 122 094707 
(2005). 

[209] J.J. Crassous et al, J.Chem.Phys. 125 204906 (2006). 
[210] M. Fuchs and M.E. Gates, J.Phys.:Condens.Matter 17, 
S1681 (2005). 

[211] N.G. van Kampen, Ghapter XV Section 3, in Stochas- 
tic Processes in Physics and Chemistry (North-Holland, 
Amsterdam, New York, Oxford 1981). 

[212] D.A. McQuarrie, Statistical mechanics (University Sci- 
ence Books, 2000). 

[213] M.E. Gates, Ann.Henri.Poincare 4, Suppl. 2 S647 
(2003). 

[214] G.H. Fredrickson and R.G. Larson, J.Ghem.Phys. 86 
1553 (1987). 

[215] M. Frank, D. Anderson, E.R. Weeks and J.F. Morris, 

J. Fluid Mech. 493 363 (2003). 
[216] H. Tresca, Gomptes Rendus Acad. Sci. Paris, 59 754 

(1854). 

[217] R. von Mises, Gottinger Nachrichten, math.-phys. 
Klasse, 582 (1913). 

[218] R. Hill The Mathematical Theory of Plasticity (Oxford 
University Press, 1971). 

[219] There is a useful short discussion of the von Mises, 
Tresca and generalized pressure dependent yield crite- 
ria in metals in section 4.2.3.3. of G.A. Schuh, T.G. 
Hufnagel and U. Ramamurty, Acta Materialia 55 4067 
(2007). 

[220] G. Ovarlez, Q. Barrel and P. Goussot, Nature Materials 
9 115 (2010). 

[221] P. Schall, D.A. Weitz and F. Spaepen, Science 318 1895 
(2007). 



